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ABSTRACT 


In  the  next  phase  of  Lincoln  Laboratory's  SWAT  (Short- Wavelength  Adaptive  Techniques) 
program,  the  performance  of  a  241 -actuator  adaptive-optics  system  will  be  measured  using  a 
variety  of  synthetic-beacon  geometries.  As  an  aid  in  this  experimental  investigation,  a  detailed  set 
of  theoretical  predictions  has  also  been  assembled  The  computational  tools  that  have  been  applied 
in  this  study  include  a  numerical  approach  in  which  Monte-Carlo  ray-trace  simulations  of 
accumulated  phase  error  are  developed  and  an  analytical  analysis  of  the  expected  system  behavior. 

This  report  describes  the  basis  of  these  two  computational  techniques  and  compares  their 
estimates  of  overall  system  performance.  Although  their  regions  of  applicability  tend  to  be 
complementary  rather  than  redundant,  good  agreement  is  usually  obtained  when  both  sets  of 
results  can  be  derived  for  the  same  engagement  scenario. 


PREFACE 


From  1982  to  1991  Lincoln  Laboratory  was  supponed  by  the  Defense  Department  to  develop 
the  technology  of  uncooperative-target  adaptive  optics.  A  principal  feature  of  this  program  was  the 
development  of  synthetic  beacons  as  the  reference  source  for  the  adaptive-optics  wavefront  sensor. 
This  technique  overcame  problems  associated  with  point-ahead  and  target  brighmess. 

In  the  late  I980’s  it  became  clear  that  synthetic-beacon  technology  could  be  of  great  value  to 
the  astronomical  community,  who  routinely  deal  with  dim  objects.  Moreover  the  basic  concept  had 
been  reinvented  by  the  astronomers  with,  apparently,  no  input  from  the  US  Defense  commumry. 
The  civilian  community  seemed  on  the  verge  of  embarking  on  a  major  development  program 
duplicating  work  the  US  Government  had  already  paid  for. 

It  seemed  incumbent  therefore  on  the  DoD  to  provide  to  the  astronomers  the  wealth  of 
information  that  had  been  developed  over  the  past  decade.  So,  in  May  1991  the  US  Air  Force 
released  virtually  all  of  the  development  work  and  associated  documentation.  The  immediate 
impact  was  that  Defense  workers  began  to  present  their  material  in  open  publication.  We  are  now 
moving  to  republish  a  number  of  documents  that  seem  of  greatest  import.  To  some  extent  these 
documents  are  dated  as  far  as  military  work  is  concerned,  but  they  arc  likely  to  be  of  great  value  to 
the  astronomers.  At  a  minimum  they  are  of  historical  and  archival  significance. 

The  report  that  follows  is  one  of  many  documents  we  are  in  the  process  of  releasing.  The 
original  report  numbered  SWP-8  (“SWAT  System  Performance  Predictions,”  dated  1  January 
1990)  is  now  updated.  If  you  possess  the  original,  it  should  be  destroyed  in  accordance  with 
current  document-control  procedures.  We  are  providing  a  copy  of  this  repor  to  the  Defense 
Technical  Information  Center  (DTIC)  so  that  users  may  request  it.  There  arc  intended  to  be  no 
distribution  limitations  on  this  document 
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1.  INTRODUCTION 


The  investigation  of  novel  beam-control  methods  for  non-cooperative  targets  has  been  a 
Lincoln  Laboratory  effon  since  1983,  when  the  tirst  phase  of  the  the  SWAT  {Shon-Wavelength 
Adaptive  Techniques)  program  was  initiated  by  the  Laboratory's  High-Energy  Laser  Beam-Contro' 
and  Propagation  Group.  The  second  phase  of  that  program  is  now  underway  at  the  Group's  60- 
cm  adaptive-optics  facility  operated  at  the  Air  Force  Maui  Optical  Station  (AMOS)  at  Maui.’  ^ 
These  experiments  have  been  generating  beam-compensation  dar  since  August,  1988. 

The  purpose  of  the  current  SWAT  experiments  is  to  validate  the  concept  of  cToaiing  synthetic 
turbulence  probes  by  exploiting  low-altitude  Rayleigh  backscatier  or  resonance  backscaiter  from 
the  sodium  layer  at  90  km.  Several  beacon  methods  have  been  proposed,  most  of  \  hich  involve 
either  the  placement  of  beams  in  the  target's  point-ahead  direction,  either  focused  beams  (the  so- 
called  "A"  method  approach)  or  laterally-sheared  collimated  beams  ("S"  method).  Lincoln 
Laboratory  has  chosen  to  concentrate  on  the  former  approach  and  to  perform  parametric  studies  of 
effects  such  as  beacon  altitude,  beacon  size,  placement  geometry,  and  reconstruction  algorithm. 

Associated  with  the  SWAT  experimental  program  is  an  active  theoretical  effort  to  simulate  and 
compare  the  performance  of  a  wide  range  of  synthetic-beacon  constructs.  These  studies  have  taken 
several  forms  but  have  recently  evolved  along  two  principal  paths  —  Monte-Carlo  ray-trace 
simulations  of  accumulated  phase  error  and  analytical  descriptions  of  expected  ensemble  behavior. 
Both  methods  are  suitable  for  the  sampling  geometries  of  immediate  interest,  and  their  domains  of 
applicability  tend  to  be  complementary.  For  those  beacon  scenarios  to  which  both  techniques  have 
been  applied,  the  resulting  predictions  arc  found  to  be  in  excellent  agreement. 

This  document  is  an  interim  report  presenting  the  more  important  aspects  of  these  on-going 
investigations.  A  brief  overview  of  the  beacon  constructs  mat  will  be  tested  by  the  SWAT  system 
is  given  in  the  next  section,  along  with  a  description  of  the  baseline  data-processing  approach.  In 
Sections  2  and  3  the  essential  characteristics  of  the  performance  models  are  outlined;  additional 
information  on  both  analysis  methods  can  be  found  in  the  Appendices. 

The  main  body  of  this  report  is  contained  in  Section  5,  where  a  series  ot  parametric 
performance  studies  are  presented.  Included  are  comparisons  as  a  function  of  source  altitude, 
number  of  beacons,  and  beacon  extent.  To  expedite  precise  correlation  with  field  data,  most  of  the 
predictions  are  given  in  terms  of  both  wavefront-control  figure-error  and  tilt-removed  far-field 
Strehl. 
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Although  all  the  present  results  are  specific  to  the  SWAT  beam-control  system,  in  a 
subsequent  report  these  analyses  will  be  generalized  to  include  apenure  dimensions  arpropnaie  for 
operanonal  high-energy  laser  systems. 


2.  FOCUSED-BEACON  SAMPLING  CONCEPTS  AND 
THEIR  IMPLEMENTATION  IN  THE  SWAT 
EXPERIMENTAL  PROGRAM 

The  earliest  proposals  incorporating  synthetic-beacons  as  a  means  of  sampling  atmospheric 
turbulence  suggested  the  use  of  a  single  pulsed  laser  to  place  a  focused  beam  at  a  nominal  altitude 
of  10  km;  at  this  range  the  Rayleigh  backscatter  from  a  5-joule  laser  is  sufficient  to  drive  a  photon- 
noise-limited  wavefront  sensor  operating  at  visible  wavelengths.  It  was  soon  apparent,  however, 
that  this  simple  construct  contains  several  fundamental  defects,  the  most  obvious  of  which  is  the 
lack  of  good  tracking  information  along  the  point-ahead  path.  A  second  concern  is  the  large 
sampling  error  that  arises  when  a  low-altitude  source  is  used  to  replace  a  physical  beacon  attached 
to  the  target.  This  source  of  error,  referred  to  as  focal  anisoplanatism,  places  a  restriction  on  the 
aperture  size  that  can  be  corrected  by  a  single  source  at  a  specified  altitude.  Figure  2-1  illustrates 
the  relevant  geometry  and  shows  the  two  primary  sources  of  anisoplanatic  error  -  the  unsampled 
atmosphere  above  the  beacon  and  the  incorrectly  sampled  atmosphere  below  the  beacon. 

Working  under  the  assumption  that  the  problems  just  described  arc  not  insurmountable,  the 
adaptive-optics  community  has  developed  a  number  of  alternative  sampling  approaches  with  a  view 
toward  obtaining  better  tracking  information  and  improving  the  figure  (tilt-removed)  compensation 
for  large  apertures.  The  SWAT  program  incorporates  an  experimental  investigation  of  both  the 
tracking  and  figure  components,  with  particular  emphasis  on  the  latter.  •  This  report  deals  strictly 
with  the  focal  anisoplanatic  effects  associated  with  single-beacon  and  multiple-beacon  deployment 
strategies. 

2.1  THE  SWAT  BEACON  GEOMETRIES 

In  Section  4.2  it  is  shown  that  the  error  variance  associated  with  focal  anisoplanatism  for  a 
single  artificial  source  can  be  expressed  by  the  relationship 

(2-1) 

where  D  is  the  aperture  diameter  and  H  is  the  beacon  altitude.  Although  one  could,  in  principle, 
achieve  an  arbitrarily  small  error  by  making  H  sufficiently  large,  in  practice  it  is  likely  to  be 
extremely  difficult  to  obtain  adequate  signal  from  a  Rayleigh  beacon  above  20  km.  An  alternative 
method  is  to  reduce  the  effective  transmitter  diameter  by  increasing  the  number  of  synthetic 
sources,  each  of  which  is  used  to  drive  an  aperture  section.  The  process  of  combining  gradient 
measurements  from  several  beacons  (often  referred  to  as  "stitching")  has  received  considerable 
attention  in  recent  years. 
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Figure  2-1.  "A"  method  beacon  geometry. 


The  SWAT  program  represents  the  community's  first  attempt  to  experimentally  verify  the 
practical  utility  of  the  synthetic-beacon  concept  and  the  feasibility  of  multiple-source  stitching  as  a 
means  of  correcting  arbitrarily  large  apertures.  The  system’s  60-cm  telescope  is  small  enough  to 
permit  a  good  correction  to  be  made  with  a  single  source,  yet  large  enough  to  show  a  measurable 
performance  difference  between  compering  multiple-beacon  methods.  The  baseline  test  sequence 
will  include  the  following  four  configurations,  each  of  which  is  studied  in  detail  later  in  this  repon. 
The  associated  beacon  geometries  are  depicted  in  Figure  2-2. 

Single  Beacon.  The  baseline  experiment  exploits  low-altitude  Rayleigh  scattering  to 
generate  a  singe  artificial  source.  Using  a  laser  capable  of  generating  a  pulse  energy 
of  the  order  of  5  Joules,  a  beacon  having  a  1-km  depth  can  be  placed  at  altitudes  as 
high  as  10  km.  The  return  radiation  is  sensed  by  a  gated  wavefront  sensor  that 
produces  gradient  data,  which  is  subsequently  processed  to  form  a  phase-front 
estimate  through  the  use  of  conventional  Gauss  reconstruction  techniques. 

Four  Beacons.  In  order  to  reduce  the  effects  of  focal  anisoplanarism,  the  aperture  is 
divided  into  four  sections  that  are  each  driven  by  a  separate  source.  In  the  SWAT 
system,  the  Rayleigh  beacons  are  placed  by  sequential  pulses  so  that  each  subapenure 
is  able  to  obtain  a  gradient  measurement  from  each  beacon.  This  arrangement  allows 
the  greatest  flexibility  in  data  processing. 

Single  Beacon  —  Dual  Altitude.  Low  spatial-frequency  anisoplanatic  errors  arising 
from  the  use  of  a  single  low-altimde  source  can  be  measured  and  corrected  in  a 
second  mirror  update  by  adding  a  high-altitude  source  that  is  measured  by  large 
subapertures.  By  combining  gradient  outputs  to  produce  four  large  subapenures  in 
this  second  iteration,  the  signal-to-noise  improvement  makes  section-tilt 
measurements  from  Rayleigh  beacons  above  15  km  feasible,  and  resonant  returns 
from  the  sodium  layer  at  90  km  can  also  be  exploited.  In  addition  to  the  six  dye 
lasers,  the  SWAT  system  may  also  incorporate  a  500-mJ  sum-frequency  Nd:YAG 
laser  for  the  purpose  of  generating  resonant  sodium  returns.^  Work  is  also 
progressing  on  a  7  J  laser  sodium  laser  that  would  provide  single-beacon  data  from 
90  km. 

Four  Beacons  —  Dual  Altitude.  In  this  scenario  multiple  low-altitude  beacons  are 
combined  with  a  single  high-altitude  source  for  optimal  correction  of  both  low  and 
high  spatial-frequency  wavefront  errors. 


Figure  2-2.  SWAT  beacon  geometries. 
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Parametric  investigations  of  such  effects  as  beacon  altitude,  radial  extent,  axial  extent,  and  later,*! 
displacement  will  also  be  included  in  this  experimental  sequence. 

The  SWAT  system  has  been  designed  to  test  a  diverse  set  of  phase  sampling  and 
compensation  methods  and  will  eventually  include  a  bank  of  six  dye  lasers  operating  at  508  nm  and 
a  Nd'.YAG  laser  capable  of  exciting  sodium-resonance  backscattcr  from  90  km.  Figure  2-3  shows 
the  major  optical  elements  of  this  system  and  illustrates  the  placement  of  the  fast  switching  mirrors 
used  to  sequentially  inject  beacons  into  the  outgoing  beam  path.  In  the  first  set  of  measurements, 
far-field  scoring  will  be  performed  by  recording  star  images  obtained  just  before  and  after  the 
application  of  a  pulsed  correction.  A  summary  description  of  the  key  wavefront  control 
components  is  prc'vided  in  Table  2-1. 

TABLE  2-1 

Wavefront-Control  System  Characteristics 

Digital  Reconstructor 

•  241  parallel  channels 

•  10^  multiply  and  add  operations/sec 

Precision  Tracking  Laser 

•  7-watt  argon-ion  laser 

•  X.=  488  nm 

Beacon  Lasers  (6) 

•  5-Joule/pulse  dye  laser 

•  X=  508  nm 

•  Beam  quality  =  20  X/D 

Beacon  Laser  (1) 

•  500-mJ/pulse  Nd;YAG 
sum-frequency  laser 

•X  =  589  nm 


Telescope 

•  60-cm  Cassegrain 

•  18-cm  central  obscuration 

Wavefront  Sensor 

•  16x16  modified  Hartmann  array 

•  256  X  and  y  gradients/sample 

•  240-iisec  readout  delay 

•  100-jisec  reconstruction  delay 

•  X  /15  performance  with  5,000  photons 

Deformable  Mirror 

•  17x17  edge-clipped  array 

•  241  active  actuators 

•  3.5-cm  acmator  spacing 

•  300-iisec  response  time 


Reference 

Star 


"■e  2-3.  Multiple-beacon  compensation  experiment. 


2.2  MULTIPLE-BEACON  PHASE  RECONSTRUCTION 

The  reconstruction  of  an  optical  wavefront  from  a  set  of  local  phase  gradient  measurements 
has  received  extensive  treatment  in  the  literature.'*  The  standard  derivation  begins  with  a  model  in 
which  the  gradient  measurements  form  a  grid  connecting  phase  points  situated  at  the  nodes  of  the 
grid.  The  raatiix  representation  of  this  model  has  the  form 


g  =  A<p  ,  (2-2) 

where  g  is  a  vector  comprising  the  x  and  y-phase  gradients,  and  (p  is  the  phase  vector.  The  A 
matrix  describes  the  physical  relationship  between  pairs  of  phase  points  and  their  connecting 
gradients.  The  inverse  of  Equation  (2-2)  provides  a  means  of  obtaining  a  best  estimate  of  the 
phase  ensemble  from  a  set  of  noisy  (and  therefore  inconsistent)  measurements  of  the  local  phase 
difference 

9>  =  {(A*A)-U‘}g  =  Bg  .  (2-3) 

It  can  be  shown  that  the  B  matrix  is  the  optimal  estimator  if  the  gradient  noise  is  random  and 
uncorrelated. 

There  are  two  major  complications  to  be  addressed  in  dealing  with  multiple  sources;  these 
arise  from  the  manner  in  which  the  beacons  are  projected  onto  the  scattering  layer  and  the  variation 
of  the  quality  of  the  gradient  measurements  as  a  funedon  of  aperture  posidon.  Both  of  these 
complicadons,  and  a  brief  description  of  the  proposed  solutions,  are  illustrated  in  Figure  2-4. 

In  order  to  project  an  array  of  sources  that  are  properly  positioned  over  their  respective 
aperture  sections,  each  laser  beam  must  be  introduced  into  the  beam  train,  tilted  by  the  correct  angle 
prior  to  transmission,  and  the  return  beam  must  be  detilted  by  the  same  angle  before  introduction 
into  the  wavefront  sensor.  Due  to  up-leg  anisoplanatic  effects,  there  is  no  beacon-projection 
method  that  will  allow  the  beam  positions  to  be  precisely  controlled.  Furthermore,  any  mistake 
made  in  the  process  of  determining  relative  source  location  will  contribute  to  a  net  "section- 
stitching"  error.  The  best  method  of  removing  this  beacon-position  sensitivity  is  to  subtract  the 
?perture-average  tilt  component  from  each  gradient  ensemble.  The  net  error  associated  with  this 
estimation  process  is  quantified  in  Section  4.3 

At  this  point  in  the  data  processing,  each  subaperture  has  one  gradient  sample  associated  with 
each  synthetic  source.  In  a  statistical  sense,  the  quality  of  these  measurements  will  be  related  to  the 
inverse  of  their  radial  distance  from  the  section  centers.  These  data  can  be  combined  through  a 
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Gradient  Error 


Problem  #1 


Solution 


Differential  Tilts  are  Introduced  Modify  the  Data  to  Remove  the 

by  Beacon  Misposition  Errors  Full- Aperture  Tilt  Component 


Problem  #2 


Solution 


Duplicate  Measurements  are  Use  Linear  Estimation  Techniques 

Available  and  Gradient  Errors  Vary  to  Develop  a  Single  Gradient  Array 

as  a  Function  of  Aperture  Position  with  Least-Square  Error  Characteristics 


Figure  2-4.  Treatment  of  the  multiple-beacon  sampling  problem. 


linear  estimation  process  defined  in  the  following  manner 

g  =  Wy  , 


(2-4) 


where  g  is  the  best  estimate  of  the  gradient  array,  y  represents  the  complete  set  of  tilt-corrected 
gradient  measurements,  and  W  is  the  optimal  weight  matrix.  The  optimal  W  matrix  can  be  derived 
from  Gauss-Markov  estimation  theory.^  The  solution  for  the  general  case  is 

W  =  (u'^c;'u)“‘u'^c;'  .  (2-5) 

where  is  the  gradient-noise  covariance  matrix  and  U  is  the  vector 

T 

1 


If  the  covariance  matrix  is  diagonal  (i.e.  the  gradient  noise  is  uncorrelated),  the  linear  weights  are 
found  to  be  inversely  proportional  to  the  gradient  variances, 
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In  the  final  computational  step,  a  best  estimate  of  the  phase  matrix,  <p,  is  obtained  from  g  through 
the  operation  described  in  Equation  (2-3). 


3.  NUMERICAL  EVALUATION  OF  SWAT  PERFORMANCE 
USING  RAY-TRACE  TECHNIQUES 

Monte  Carlo  simulations  incorporating  ray-trace  techniques  have  frequently  been  used  in  the 
study  of  turbulence-induced  phase  distortion  processes.  Although  this  method  has  been  criticized 
as  being  inelegant,  it  can  accommodate  virtually  any  effect  that  is  relevant  to  the  performance  of  an 
adaptive  phase-correction  system.  In  addition,  estimates  of  Strehl  can  be  obtained  directly  from 
far-field  beam  profiles  derived  from  the  Fourier  transform  of  the  near-field  phase. 

The  accuracy  of  the  output  generated  by  a  numerical  simulation  is  strongly  dependent  upon  the 
programmer's  ability  to  properly  incorporate  all  important  physical  effects,  to  create  a  model  with 
sufficient  sampling  resolution,  and  to  establish  a  large  enough  ensemble  to  provide  statistically 
meaningful  results.  All  these  requirements  must  be  weighed  against  the  attendant  costs  of 
computational  memory  and  execution  time.  The  programming  approach  described  in  this  section 
was  structured  to  satisfy  all  these  constraints. 

3.1  NUMERICAL  CODE  STRUCTURE 

3.1.1  Source  Construction  and  Turbulence  Profile 

The  structure  of  the  propagation  code  ELVIRA  used  in  this  study  of  the  SWAT  experimental 
system  is  outlined  in  Figure  3-1,  The  logical  flow  begins  with  an  atmospheric  model  represented 
by  a  series  of  20  phase  screens,  appropriately  weighted  and  positioned  between  ground  level  and 
20  km.  In  order  to  bracket  the  expected  performance  over  a  reasonably  wide  range  of  turbulence 
conditions,  both  the  HV-21  and  SLC-Day  models  are  incorporated  in  the  code;  these  turbulence 
profiles  are  defined  in  Appendix  E,  and  the  associated  screen  weights  and  altitudes  are  listed  in 
Table  3-1. 

Each  of  the  screens  is  represented  by  a  high-rcsolution  phase-distortion  array  covering  a  76  x 
76  grid  with  0.875-cm  spacing.  This  allows  the  subapertures,  which  are  assumed  to  have  a  3.5- 
cm  spacing,  to  be  represented  by  a  4  x  4  subarray.  It  is  the  incorporation  of  this  fine  structure  in 
the  simulation  that  permits  gradient-sampling  and  influence-function  effects  to  be  investigated.  Ten 
turbulence  realizations  are  applied  to  each  sampling  scenario;  the  associated  coherence  length  (r^) 
and  isoplanatic  angle  (&„)  parameters  are  listed  in  Table  3-2.  In  all  cases  the  ensemble  average  is 
within  5%  of  the  expected  value  for  the  specified  turbulence  profile. 
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TABLE  3-1 

Phase  Screen  Altitudes  and  Weights  for  Ray-Trace  Simulation 


Screen  Altitude  (km) 

HV-21  Phase  Screen 
Weights  ) 

SLC-Day  Phase  Screen 
Weights  ) 

0.125 

1.623x10*''^ 

8.236x10*’^ 

0.375 

1.807x10'^^ 

3.250x10*’^ 

0.625 

5.506x10'’^ 

3.250x10*’® 

0.875 

3.857x10'^'* 

2.838x10*’® 

1.250 

5.902x10*’"^ 

2.245x10*’® 

1.750 

4.224x1 0'’"^ 

7.954x10*’^ 

2.250 

3.027x1 0'^^ 

3.714x10*’"’ 

2.750 

2.171x10'^^ 

2.032x10*’"’ 

3.250 

1.563x10’’^ 

1.232x10*’"* 

3.750 

1.137x10'^^ 

8.037x10*’® 

4.250 

8.485x10*’® 

5.534x10*’® 

4.750 

6.648x10*’® 

3.973x10*’® 

5.250 

5.629x10*’® 

2.950x10*’® 

5.75C 

5.240x10*’® 

2.251x10*’® 

6.250 

5.317x10*’® 

1.757x10*’® 

6.750 

5.712x10*’® 

1.397x10*’® 

7.250 

6.288x10*’® 

1.223x10*’® 

7.750 

6.926x10*’® 

1.173x10*’® 

9.000 

3.235x10*’"’ 

4.284x10*’® 

15.000 

7.647x10*’"’ 

1.586x10*’"* 

TABLE  3-2 

Turbulence  Parameters  for  Simulation  Test  Ensemble 


HV-21 

Model 

SLC-Day  Model 

Run  Number 

fo  (cm) 

do(lirad) 

To  (cm) 

eo(lirad) 

1 

4.2 

7.8 

4.4 

2 

6.3 

6.5 

5.8 

■H 

3 

4.8 

7.0 

4.9 

12.5 

4 

4.6 

6.2 

5.2 

12.7 

5 

4.4 

5.2 

4.8 

10.8 

6 

5.1 

5.2 

5.2 

11.1 

7 

4.6 

6.6 

4.6 

12.5 

8 

4.8 

8.8 

5.1 

14.9 

9 

6.1 

9.4 

5.4 

15.2 

10 

5.3 

5.3 

5.7 

11.1 

Average 

5.0 

6.8 

5.1 

12.6 
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3.1.2  Ray-Trace  Procedure 

As  indicated  in  Figure  3-2,  two  sets  of  ray-trace  calculations  are  performed  through  the  phase- 
screen  ensemble.  The  reference  or  "truth”  wavefront  is  generated  by  summing  the  array  points 
along  collimated  ray  paths.  In  a  similar  fashion,  the  wavefronts  associated  with  the  low-altitude 
sources  are  obtained  by  adding  the  encountered  phase  points  {with  appropriate  interpolation)  along 
the  slanted  rays  connecting  the  source  and  the  measurement  point  at  the  collection  apenure. 
Separate  high-resolution  wavefront  fields  are  produced  for  each  of  the  synthetic  beacons  used  in  a 
given  sampling  scenario. 

3.1.3  Derivation  of  Phase  Gradients 

Once  ail  of  the  phase  arrays  have  been  defined  in  the  plane  of  the  receiver  apenure.  a 
measurement  model  can  be  applied  to  extract  simulated  gradient  data.  The  SWAT  sensor 
incorporates  a  "toe-to-tail"  sampling  arrangement  that  requires  separate  x  and  y  cameras.  The 
spatial  relationship  between  the  actuator  position  and  the  gradient  sampling  regions  is  described  in 
Figure  3-3.  Individual  subaperture  tilts  over  the  enclosed  4x4  subarrays  are  determined  by  fitting 
a  least-square  plane  to  the  included  phase  points.  The  gradient  output  array  is  formed  from  those 
measurements  that  fall  within  the  60-cm  SWAT  aperture. 

3.1.4  Phase  Reconstruction 

The  general  approach  for  obtaining  a  best  estimate  of  the  wavefront  error  was  presented  earlier 
in  Section  2,2.  Both  the  ray-trace  and  analytical  analyses  described  in  this  rcpon  use  a  somewhat 
simplified  approach  for  combining  data  from  multiple  synthetic  sources,  in  which  the  final  gradient 
matrix  is  constructed  from  measurements  generated  by  the  nearest  beacon.  Specifically,  the 
reconstructed  phase  is  represented  by  the  matrix  operation 

<p  =  BSTg  .  (3-1) 

where  g  is  a  gradient  vector  made  up  of  measurements  from  all  available  sources,  T  is  the  matrix 
that  removes  overall  tilt  from  each  of  the  sections,  and  S  is  the  operation  of  selecting  the  most 
accurate  gradient  sample  from  the  available  data.  (The  B  matrix  is  the  same  one  used  previously 
for  single-source  reconstruction.)  Direct  comparisons  of  single  and  multiple-beacon  configurations 
can  be  performed  on  the  same  set  of  turbulence  realizations,  and  competing  stitching  aigonthms 
can  be  readily  evaluated.  At  this  point  in  the  analysis,  an  estimate  of  the  anisoplanatic  error 
associated  with  the  beacon  geometry  can  be  made  by  comparing  the  collimated  phase  at  the  actuator 
points  with  the  reconstructed  phase  array. 
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Phase  Screens 
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Figure  3-2.  Ray-trace  geometry. 


Sample  Point  Actuator  Position  Gradient  Sample 
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Figure  3-3.  Gradient  sampling  of  the  high-resolution  grid. 


3.1.5  Deformable  Mirror  Influence  Function  Effects 

The  final  step  in  the  simulation  of  phase  compensation  occurs  when  an  influence  function  is 
applied  to  the  actuator-drive  positions  to  create  a  high-resolution  figure  over  the  mirror  surface. 
This  process  can  be  represented  as  a  convolution  process  having  the  form 

M{i,j)  =  F{i,j)®5{i-An,j-Am)  ,  (3-2) 

where  F(i,  j)  is  the  mirror  influence  function,  and  S(i-4n,  j-4m)  places  a  delta  function  at 
every  actuator  position.  For  convenience,  the  actuators  are  assumed  to  be  infinitely  rigid  so  that 
F(/,  j)  is  zero  when  {i  =  4n,n^  0}  and  {7  =  4m,  m^O}. 

Figure  3-4  illustrates  the  two  influence  functions  that  have  been  used  in  this  study,  each  of 
which  has  the  characteristic  of  correctly  replicating  pure  tilt  and  piston.  The  function  labeled 
"SWAT  mirror”  has  been  derived  from  measurements  of  the  faceplate's  response  to  the  motion  of  a 
single  actuator,  but  the  attendant  fitting  error  differs  little  from  the  standard  trapezoidal  form 
defined  below 


F{-  •) -  l^~'^  J  ~M  :  for  I  /  i  <  4  and  ! 7  1  <  4  ,  and 
16 

=  0  :  elsewhere  .  (3-3) 

For  simplicity,  the  trapezoidal  function  has  been  applied  in  most  of  the  cases  described  in  this 
report.  The  standard  simulation  also  ignores  the  central  obscuration  present  in  the  SWAT  mirror 
and  derives  the  phase  error  function  from  just  the  inner  55  cm  of  the  mirror's  diameter  in  order  to 
eliminate  anomalous  edge  effects. 

3.1.6  Performance  Evaluation 

Once  the  mirror  surface  figure  has  been  generated  in  the  manner  just  described,  an  error 
function  representing  the  difference  between  the  collimated-source  phase  and  the  mirror  phase  is 
established.  The  overall  piston  term  is  then  removed  from  this  result,  and  the  apenure  tilt  and 
figure  components  are  separated  by  removing  overall  tilt  from  the  difference  function.  For  figure 
values  smaller  than  about  1.5  rad^,  a  reasonably  accurate  estimate  of  the  tilt-removed  Strehl  can  be 
obtained  from  the  extended  Marechal  approximation 

SrreW  =  exp(-<T^)  .  (3-4) 

Alternatively,  the  Strehl  ratio  can  be  computed  directly  from  the  far-field  image  derived  from  the 
Fourier  transform  of  the  figure  error  function.  To  avoid  aliasing  errors  the  code  embeds  the 
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Trapezoidal  Influence  Function 


Best  Estimate  of  the  SWAT  n  Mirror  Influence  Function 


Figure  3-4.  Irrfluencefiincnons  used  in  the  numerical  ray-trace  simulation. 
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apenure  function  in  a  512  x  512  field  of  zeros  prior  to  computing  the  Fourier  transform.  Figure 
3—5  shows  an  example  of  a  near-field  phase  error  function  and  its  associated  far-fieid  image 
computed  in  this  manner.  Obtaining  accurate  Strehl  ratios  for  highly  scintillated  beams  is  a  non¬ 
trivial  exercise;  the  process  incorporated  both  in  this  study  and  the  analysis  of  field  data  is 
described  in  Appendix  D.  Final  performance  values  are  quoted  on  the  basis  of  average  phase 
variance  or  average  Strehl  over  the  sample  ensemble  of  10  turbulence  realizations. 

3.2  INTRODUCTORY  SYSTEM-PERFORMANCE  OBSERVATIONS 

A  number  of  general  statements  can  be  made  regarding  the  interaction  of  cenain  types  of 
system-induced  errors  and  the  relationship  between  tilt-removed  error  variance  and  Strehl.  These 
observations  appear  to  be  valid  for  a  wide  range  of  data-coUection  scenarios. 

3.2.1  Gradient  Sampling  and  Mirror  Fitting  Effects 

The  numerical  code  has  been  designed  to  include  a  number  of  important  hardware-related 
effects  that  are  not  easily  modeled  analytically.  In  particular,  the  spatial  sampling  characteristics  of 
the  wavefront  sensor  and  the  fitting  error  introduced  by  the  deformable  mirror  are  explicitly 
included  in  the  model. 

Reconstruction  techniques  that  produce  phase  estimates  from  gradient  data  usually  include  the 
implicit  assumption  that  each  measured  gradient  is  equivalent  to  a  phase  difference  between 
adjacent  actuator  locations.  As  shown  earlier,  in  Figure  3-3,  this  assumption  is  only  approximately 
valid,  since  each  measurement  is  made  over  a  finite  spatial  region.  By  comparing  the  reference 
wavefront  sampled  at  actuator  grid  points  with  the  reconstructed  phase  values  for  an  ensemble  of 
turbulence  realizations,  it  was  determined  that  the  average  variance  associated  with  the  sampling 
geometry  is  0.09  rad^.  This  value  is  small  compared  to  most  other  sources  of  error  encountered  in 
these  simulations  and  also  appears  to  be  highly  correlated  with  other  high  spatial-frequency  effects. 

Fitting  error  refers  to  a  deformable  mirror's  ability  to  correctly  reproduce  the  high  spatial 
frequencies  associated  with  Kolmogorov  turbulence.  For  an  aperttne  of  infinite  extent,  the  fitting 
error  is  predicted  to  have  the  approximate  form^ 

a%,^l  =  0.32(<i,^  / r,f'\  (3-5) 

where  is  the  subaperture  diameter.  For  the  nominal  SWAT  test  conditions  cm  and 
(^sub  -  3-5  cm,  which  yields  a  fitting  error  prediction  of  0. 18  rad^. 


Near-Field  Figure  Error 


Far-Field  Strehl 


Figure  3-5.  Numerical  simulation  sample  output. 


The  SWAT  deformable  mirror  contains  a  total  of  341  active  acruatorr,  of  which  241  fall  within 
the  nominal  60*cm  aperture  dianxiter.  The  remaining  actuators  form  a  guard  ring  that  can  be  held  at 
the  average  drive  voltage  as  a  means  of  isolating  the  active  section  from  the  face-sheet  support 
structure.  Using  the  HV-21  turbulence  model  to  generate  an  ensemble  of  input  wavefronts,  fitting 
error  effects  were  measured  over  concentric  mirror  regions  between  40  and  60  cm.  As  shown  in 
Figure  3-6,  the  tilt-removed  variances  for  both  the  trapezoidal  and  SWAT  influence  functions  are  in 
accordance  with  the  predicted  value  out  to  55  cm. 

3.2.2  Comparison  of  Anisoplanatic  Error  and  Total  Near-Field  Error 
Variance 

The  validity  of  error  estimates  obtained  by  summing  variances  from  multiple  error  sources  has 
frequently  been  questioned  in  adaptive-optics  problems,  particularly  when  the  noise  mechanism  is 
associated  with  unconected  high-spatial-frequency  turbulence  components.  This  situation  is 
illustrated  in  Figure  3-7,  which  compares  the  total  phase  error  predicted  by  the  code  for  a  single- 
beacon  geometry  at  a  range  of  altitudes  with  the  error  from  focal  anisoplanatic  effects.  At  low 
altitudes,  where  focal  anisoplanatic  effects  are  most  significant,  the  two  values  are  nearly  identical. 
The  fitting  error  component  of  0.2  rad^  only  appears  when  the  anisoplanatic  error  approaches  zero. 
Over  the  region  of  greatest  interest,  where  the  error  variance  is  between  0.5  and  1.5  rad^,  system 
effects  generally  add  no  more  than  0. 1  rad^  (k/20)  to  the  anisoplanatic  error.  This  observation  will 
make  it  easier  to  relate  the  analytical  results  described  in  the  next  section  to  overall  system 
performance 

3.2.3  Comparison  of  Error  Variance  and  Far-Field  Strehl 

The  normal  treatment  of  system  performance  will  present  data  in  terms  of  error  variance.  One 
assumes  that  this  information  can  be  translated  in  some  quantity  that  relates  to  the  far-field  beam, 
and  the  extended  Marechal  approximation  is  generally  applied  to  provide  this  comparison.  For  the 
results  of  interest  in  this  repon,  the  values  arc  most  often  quoted  for  tilt-removed  calculations,  and 
the  applicability  of  the  extended  Marechal  approximation  for  this  case  is  not  immediately  obvious. 

It  has  been  understood  for  some  time  that  small  phase  fluctuations  in  the  apenure  plane  due  to 
poor  phase  compensation  produce  a  far-field  beam  with  a  nearly  diffraction-limited  central  lobe  but 
with  a  reduced  peak  intensity.  The  lost  energy  is  typically  due  to  high  spatial  frequencies  that 
broadly  distribute  energy  in  the  far  field.  As  the  phase  is  further  corrupted,  the  beam  profile 
becomes  more  irregular,  making  a  determination  of  Strehl  much  more  difficult.  Appendix  D 
describes  an  analysis  technique  that  was  developed  in  the  course  of  studying  SWAT  data  that 
makes  extensive  use  of  histogram  transformations  to  perform  background  normalization  and 
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FITTING  ERROR  (rad 


Figure  3-6.  Fitting  error  results  as  a  function  of  aperture  mask  diameter.  The  edge  of  the  60-cm 
mask  lies  just  beyond  the  last  ring  of  active  actuators. 
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Figure  3-7.  Comparison  of  numerical  results  for  focal  anisoplanatic  error  and  total  system 

error. 
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specify  Strehl.  To  a  good  approximation,  the  Strehl  value  obtained  in  this  manner  is  the  ratio  of 
the  peak  pixel  intensity  to  the  expected  peak  intensity  for  a  diffraction-limited  beam  having  the 
same  total  energy. 


ELVIRA  uses  this  same  approach  in  obtaining  Strehl  from  the  simulated  far-field  data.  The 
algorithm  removes  overall  tilt  by  defmition  and  therefore  is  appropriate  for  this  study.  Figure  3-8 
proddes  a  comparison  of  Strehl  with  the  near-field  error  variance  for  two  types  of  noise,  the 
circles  relate  to  focal  anisoplanadc  error  introduced  by  placing  a  single  beacon  at  a  variety  of 
altitudes;  the  crosses  represent  randor  Gaussian  phase  error.  The  extended  Marechal 
approximation  is  seen  to  be  valid  for  b-  '  types  of  data  for  figure  errors  less  than  1.5  rad^ 
(corresponding  to  a  Strehl  of  about  0.22).  For  larger  phase  distonions  the  far-field  performance 
will  be  underestimated  by  the  Marechal  formula,  and  a  better  fit  to  the  code  output  can  be  obtained 
firom  the  empirical  expression 


Strehl  - 


exp^-CT^) 

O.5exp^-O.538cr^ ) 


<1.5 

:£r^^l.5 


(3-6) 


which  yields  reasonable  estimates  for  figure  errors  as  large  as  5  rad^. 


3.2.4  Correlation  of  Strehl  With  and  0^ 

Figure  erroR  from  focal  anisoplanadc  effects  are  not  simple  functions  of  low-order  turbulence 
moments  but  are,  instead,  expressed  as  infinite  Taylor  series  of  fractional  moments  (moments 
calculated  from  ground  level  to  the  beacon  altitude).  In  analyzing  future  field  data,  it  will  be 
important  to  undeRtand  the  puise-to-pulse  performance  variation  that  is  to  be  expected  for  a  given 
atmospheric  model  and  to  establish  the  correlation  (if  any)  between  compensation  performance  and 
measurable  turbulence  parametcR. 


Figure  3-9  plots  Strehl  as  a  function  of  and  6^,  for  ten  HV-21  and  SLC-Day  profile 
realizations  using  a  5-km  single-beacon  geometry.  Both  sets  of  realizations  show  a  scatter  of  about 
a  factor  of  2  in  the  estimated  correct  Strehl  and  yield  correlations  that  are  model  specific  and 
relatively  weak.  (When  the  two  turbulence  profiles  ensembles  were  separately  processed,  a 
correlation  coefficient  of  about  0.5  was  obtained  for  each  parameter.)  Note  that  the  expected  6^ 
means  of  7  and  12  firad  were  obtained  for  the  HV-21  and  SLC-Day  models  respectively,  but  that 
each  produced  the  same  average  Strehl. 
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Figure  3-8.  Comparison  offar-field  Strehl  measurement  with  the  extended  Marechal  approxi¬ 
mation  applied  to  the  near-field  phase  error.  The  circles  indicate  a  range  of  source  placement 
scenarios  that  introduce  focal  anisoplanatic  error,  and  the  crosses  relate  to  calculations  in 
which  random  phase  errors  have  been  imposed. 


28 


STREHL  STREHI 


0  S  10  1S  20 

00  ((irad) 

Figure  3-9.  Correlation  ofStrehl  estimate  with  r^  and  6^  values  for  10  turbulence  realizations 
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This  simulation  shows  that  overall  performance  is  not  well  correlated  with  either  or  6^,  and 
that  neither  of  these  parameters  is  likely  to  be  a  good  predictor  of  the  achievable  far-field  Strehl  for 
an  individual  synthetic-beacon  experiment.  This  is  not  to  imply  that  these  quantities  are 
unimportant  in  the  subsequent  data  analysis,  but  rather  that  the  error  sources  associated  with  finite- 
altirade  beacon  arrays  are  too  complex  to  be  characterized  by  a  single  turbulence  moment.  Methods 
for  obtaining  adequate  turbulence-profile  information  are  actively  sought,  but  in  the  interim  the 
performance  curves  given  in  this  repon  are  likely  to  serve  as  the  most  reliable  guide  for  the  SWAT 
experiments. 
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4.  ANALYTIC  EVALUATION  OF  SWAT  PERFORMANCE 
USING  MELLIN  TRANSFORM  TECHNIQUES 

In  Section  3  of  this  report,  a  method  was  described  that  has  been  used  to  evaluate  the 
performance  of  the  SWAT  system  through  the  application  of  conventional  ray-trace  techniques. 
The  analysis  of  much  larger  systems  can  present  severe  computational  difficulties,  however,  and 
for  this  reason  useful  scaling  laws  for  aperture  diameter  effects  are  not  easily  derived  in  this 
manner. 

In  this  section  a  different  and  more  general  approach  to  the  system-analysis  problem  is 
presented.  This  method  is  computationaily  efficient  and,  for  small  systems,  has  provided  results 
that  are  close  to  those  obtained  with  the  ray-trace  code.  The  technique  can  be  extended  to  larger 
systems  with  no  difficulty,  and  because  the  technique  breaks  the  problem  down  into  vanous 
components  one  can  readily  identify  the  major  sources  of  system  degradation. 

4.1  COMPUTATIONAL  MODEL 

The  model  used  to  calculate  the  phase  variance  assumes  the  variance  is  composed  of  two 
major  parts,  as  shown  in  Figure  4-1.  The  first  part  is  the  focal  anisoplanatic  error,  which  of  itself 
has  the  phase  variance  from  the  unsensed  turbulence  above  the  beacons  plus  the  error  made  in 
sensing  the  turbulence  below  the  beacons.  If  there  were  only  one  beacon,  the  sum  of  these  two 
erron  with  piston  and  tilt  removed  would  give  the  total  figure  variance.  For  the  case  of  more  than 
one  beacon  there  is  a  tilt  error  made  in  measuring  the  turbulence  below  the  beacons;  this  error  can 
be  different  for  each  section.  In  the  stitching  operation  of  putting  together  the  phases  of  each 
section  these  tilt  errors  combine  to  increase  the  total  error.  The  tilt  error  is  composed  of  two  parts; 
the  first  is  from  the  tilt  difference  between  the  collimated  beam  and  the  focused  beacon  beam;  the 
second  is  from  the  misplacement  of  the  beacons  above  the  center  of  each  section.  It  will  be 
assumed  that  the  relative  beacon  positions  are  measured  with  the  full  apenure  of  the  receiver. 
Because  the  ray  paths  through  the  atmosphere  are  different  for  the  various  beacons  there  will  be  a 
measurement  error  associated  with  the  evaluation  of  the  relative  beacon  positions.  To  obtain  the 
composite  phase  profile,  the  tilt  jitters  for  all  sections  are  stitched  together.  The  correlation  of  the 
tilts  affects  the  phase  variance  averaged  over  the  entire  apenure  when  the  full-apenure  tilt  is 
removed. 

Zemike  components  o^her  than  tilt  are  also  correlated  between  sections.  This  correlation  is 
neglected  with  little  resultant  error,  since  these  terms  represent  less  than  15%  of  the  total  variance 
and  fall  off  faster  than  tilt  as  a  function  of  separation. 
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Figure  4-1 .  Contributions  to  the  phase  varian.  e. 


Each  of  the  components  mentioned  above  will  be  evaluated  separately  by  using  Mellin 
transform  techniques  combined  with  the  use  of  the  appropriate  filter  functions.  Analytic 
expressions  will  first  be  found  for  the  components  of  the  focal  anisoplanatism  with  piston  and  tilt 
removed  above  and  below  the  beacon  altitude.  The  tilt  component  below  the  beacon  will  also  be 
found  and  used  in  the  stixhing  analysis.  It  will  be  assumed  in  this  analysis  that  the  sections  are 
circular  in  shape  in  order  to  get  analytic  expressions  for  the  results.  In  actuality  the  butted  sections 
are  square.  The  error  resulting  from  this  assumption  should  be  small. 

To  stitch  the  individual  sections  together  a  least-squares  estimation  procedure  is  used.  This 
approach  is  similar  to  that  of  stitching  together  the  tilt  errors  measured  in  a  wavefront  sensor  in  an 
adaptive-optics  system  to  obtain  the  phase  profile.**  The  difference  in  this  case  is  that  the  tilt  errors 
can  be  correlated;  this  modifies  the  formula  for  the  error  propagator.  Because  the  correlation  of  the 
tilts  can  greatly  affect  the  answer  it  is  necessary  to  derive  expressions  for  the  correlation 
coefficients  of  the  focal  anisoplanatic  tilts.  From  these  correlation  coefficients  the  error  propagator 
for  stitching  of  the  tilts  is  found.  The  error  propagator  multiplied  by  the  tilt  phase  variance  is  the 


stitching  component  of  the  phase  variance.  It  is  added  to  the  tilt-removed  phase  variance  of  a 
single  section  to  obtain  the  total  phase  variance. 

Finally,  the  phase  variance  from  focal  anisoplanatism  and  stitching  are  added  together  to  get 
the  total  phase  variance.  This  computation  is  performed  for  various  altitudes. 


4.2  FOCAL  ANISOPLANATISM 

In  a  previous  report  (TR-807)^  a  technique  was  developed  in  which  the  phase  variance  could 
be  written  down  by  inserting  a  filter  function  appropriate  to  a  given  problem  into  a  general 
expression  for  the  phase  variance.  Assuming  that  inner  and  outer  scale  effects  of  turbulence  are 
negligible,  the  expression  for  the  phase  variance  from  TR-807  is 


L 

c^  =  Q.2mki^dzCl{z)\dKK~^^'^F{K,z)  ,  (4-1) 

0 

where  C^(z)  is  the  mrbulence  strength,  =  2;r/ A ,  A  is  the  wavelength  of  the  propagating 
beam,  and  L  is  the  range  to  the  beacon.  The  filter  function  for  the  total  focal  anisoplanatism  below 
the  beacon  is  given  in  Equation  (3.4.49)  of  TR-807  as 


F(k,z)  =  2 


\  2L 
kD;Z 


2L 


(4-2) 


where  Dg  is  the  section  diameter. 


The  filter  function  for  the  tilt  is  the  square  of  the  difference  between  the  complex,  tilt  filter 
functions  of  a  collimated  and  focused  beam,  which  is 


F{k,  z)  = 


KDgl2 


AJ2{KDg(\-zlL)l2) 

KDgi\-zlL)t2 


(4-3) 


In  Appendix  A  the  variance  due  to  these  components  is  found.  If  the  beacon  is  below  the  top  of 
the  turbulence,  the  partial  moment  of  the  turbulence  must  be  used  instead  of  the  full  moment.  The 
fuU  moment  is  defined  as 


/i„  =  sec"^‘(^)JdzQ^(z)z”  . 
0 


(4-4) 
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Notice  that  the  zenith  angle  dependence  is  contained  in  the  uubulence  moment;  many  authors  do 
not  include  the  zenith  dependence  in  this  definition.  The  upper  moment  is  defined  as 


/tJ(L)  =  sec"^‘(.^)jcfzC^(z)z'’  , 


H 


where  H  =  Lj scc{^)  is  the  beacon  altitude.  The  lower  moment  is  defined  as 


H 


=  jdzCi{z)z'' 


(4-5) 


(4-6) 


The  effect  of  zenith  angle  is  also  incorporated  into  the  partial  turbulence  moments. 

One  finds  that  the  total  phase  variance  due  to  the  turbulence  below  the  beacon  is 

v5/3 


of  - 


The  variance  due  to  tilt  is 


0.355  -0.484  0.543 


l}  If  L*  If 


The  piston  variance  is 


oi  = 


(4-7) 


(4-8) 


0. 0835  0. 0 1 378  -i-  0. 00697 -  0. 059 1  0. 045 


(4-9) 


The  variance  with  piston  and  tilt  removed  is 


-0.50^3*2 


_  0.877  +  0.94  -  2.014^^^ 

r5/3  j2  fi  r4 


1. 204  -  0. 322  ^^4^ 


(4-10) 
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For  changes  in  the  zenith  angle,  if  the  beacon  is  kept  at  the  same  altitude,  it  is  clear  from  the 
variation  of  each  of  the  terms  in  brackets  that  the  zenith  dependence  of  the  variance  is  simply 
sec(i^).  To  find  the  effect  of  the  turbulence  above  the  beacon  one  must  use  the  filter  function  that 
removes  piston  and  tilt  which  is 


F{k,z)  =  \- 


■2y,{x:D/2)1 

2 

'4J2{KD/2)f 

kD/2 

kDI2 

(4-11) 


Using  Mellin  transform  techniques  one  can  find  the  piston  and  tilt  removed  variance  due  to  the 
turbulence  above  the  beacons  as 


=  0.0569ifej£)5/Vo(^)  •  (4-12) 

In  this  formula  the  total  aperture  diameter  must  be  inserted.  The  contribution  from  the  unsensed 
turbulence  increases  with  increasing  aperture  size,  and  one  finds  that  the  beacon  altitude  must  be 
increased  for  the  larger  diameters  in  order  to  keep  this  component  manageable. 

In  Figiue  4-2  the  contributions  ftom  the  turbulence  below  and  above  the  beacon  and  their  total 
are  plotted  for  a  single  beacon  operating  with  a  60-cm  diameter  aperture  with  the  HV-21  turbulence 
model.  The  total  variance  is  the  amount  that  one  expects  from  a  single  beacon  system. 

In  Figure  4-3  the  same  results  are  plotted  for  the  four  beacon  case.  In  this  case  the  section 
diameter  has  been  reduced  by  a  factor  of  2,  which  causes  the  lower  contribution  to  the  variance  to 
be  decreased  by  a  factor  of  3. 17. 

The  upper  contribution  is  the  same  in  the  2  cases.  To  find  the  total  variance  requires  the 
addition  of  the  stitching  error  to  the  focal  anisoplanatic  contribution. 
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Figure  4-2.  Predicted  SWAT  phase  variance  as  a  Junction  of  beacon  altitude  for  a  single¬ 
beacon  system  with  HV-21  turbulence.  The  turbulence  components  below  and  above  the 
beacon  are  given,  along  with  the  total  phase  variance. 
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Figure  4-3.  Predicted  SWAT  phase  variance  over  a  single  section  for  a  system  with  four 
synthetic  beacons  and  HV-21  Turbulence.  The  turbulence  components  below  and  above  the 
beacon  are  given,  along  with  the  total  phase  variance.  Stitching  effects  are  not  included. 
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4.3  BEACON  POSITION  MEASUREMENT  ERROR 

If  more  than  one  beacon  is  used,  the  tilts  of  the  individual  sections  must  be  measured.  Errors 
in  the  measurements  are  important  because  the  tilt  component  of  turbulence  is  87%  of  all  the 
turbulence -induced  phase  variance  with  piston  excluded.  There  are  two  errors  in  measuring  the 
tilt  The  first  results  from  the  difference  in  ray  paths  between  a  collimated  and  focused  beam.  This 
error  was  calculated  in  the  last  section.  The  other  result  from  the  beacons  not  being  directly  over 
the  centers  of  the  sections.  Here  it  will  be  assumed  that  the  beacons  are  projected  up  into  the  sky 
by  some  system  and  that  their  relative  positions  are  measured  by  the  full  aperture  of  the  system. 
Because  of  the  difference  in  paths  between  the  measurement  rays  through  the  turbulence,  there  will 
be  a  jitter  in  the  apparent  relative  position  of  the  beacons.  This  problem  is  analyzed  in  Appendix  B. 
The  filter  function  to  find  the  tilt  difference  between  focused  displaced  rays  that  are  separated  in 
space  is  the  product  of  the  tilt  filter  function  and  the  one  for  anisoplanatism  given  by 


F{ic,z)  = 


16 

koD 


^72(>cD[1-z/L]/2)\ 

k'D[1-2/L]/2  J 


2  l-COs|^  <5(2)11  . 


(4-13) 


Using  this  filter  function  the  tilt  variance  is  found  to  be 


12.16  ^ 
^51/3 

*,/2 
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(4-14) 


The  transition  altitude  between  the  low  and  high  altitude  solutions  is  given  by 

L 

^'"6/D-t-l  ■ 


(4-15) 


In  Figure  4-4  both  the  variance  from  this  tilt  and  the  tilt  component  of  focal  anisoplanatism 
are  plotted.  Tlie  focal-anisoplanatic  tilt  is  considerably  larger  than  that  due  to  misplacement. 
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4.4  CORRELATION  FUNCTION  OF  THE  TILT  COMPONENT  OF  FOCAL 
ANISOPLANATISM 

In  calculating  the  error  in  stitching  together  the  tilt  errors  of  the  individual  sections,  the 
correlation  of  the  tilt  is  very  important  in  determining  the  resulting  phase  variance  with  the  tilt 
removed.  For  instance,  in  Figure  4-5  are  shown  the  results  from  stitching  together  perfectly 
correlated  tilts,  and  tilts  that  arc  partially  correlated.  In  the  perfectly  correlated  case,  the  resultant 
phase  is  composed  solely  of  tilt,  and  the  tilt-removed  variance  is  zero. 


CORRELATED  PARTIALLY  CORRELATED 


FIGURE  ERROR 


Figure  4-5.  Stitching  difference  between  correlated  and  uncorrelated  tilts. 


For  turbulence  the  correlation  function  of  parallel  components  of  tilt  is  a  function  of  the  angle 
of  displacement  with  respect  to  the  tilt.  This  correlation  function  as  the  angle  varies  is  calculated  in 
Appendix  C  and  plotted  in  Figure  4-6  for  parallel,  perpendicular,  and  45  degree  displacements. 
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Figure  4-8.  Aperture  model  for  section  stitching  for  a  four-beacon  system. 


It  is  assumed  that  the  signal-to-noise  ratio  is  infinite,  and  the  measurement  noise  is  due  to  the 
phase-variance  error  incurred  in  propagating  through  turbulence.  Each  gradient  is  the  average  of 
the  difference  of  the  four  phase  points  around  that  gradient,  and  is  given  by 
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The  noise  correlation  matrix  of  the  gradients  is 

C„=(n^nJ),  (4-18) 

and  it  is  explicitly  equal  to 

1  i  p  q  0  0  0  r 

i  I  q  p  0  0  r  0 

p  q  I  i  0  r  0  0 

q  p  i  1  r  0  0  0 

Q  Q  Q  r  \  i  p  q 

0  0  r  0  /  \  q  p 

0  r  0  0  p  1  i 

r  0  0  0  q  p  i  1 

where  i=  0.558,  p  =  0.778,  q  =  0.594,  and  r  =  0.109.  These  values  come  from  the  tilt 
correlation  curves  given  in  Figures  4-6  and  4-7. 

The  least  squares  solution  for  the  phase  is 

0  =  (a^C;^ a)^ A^C;'m  =  Lm  .  (4-20) 

where  the  superscript  -t-  is  the  symbol  for  the  generalized  inverse  that  is  necessarv'  because  .he 
matrix  has  some  singular  values. 


(4-17) 


The  tilt  and  piston  removed  phase  is 


0,pf(  =  A  (a^C;^ a)^  A^C^’m  =  A  Lm  .  (4-2 1 ) 


L  is  the  tilt  removal  matrix  that  is  a  combination  of  x  and  y-tilt  removal  given  by 


(4-22) 

where 

=0.5[-l  0  1-10  1-10  1], 

(4-23) 

and 

/y  =0.5[  -1  -1  -1  0  0  0  111]. 

(4-24) 

The  piston  does  not  have  to  be  explicitly  removed  since  the  A  matrix  already  has  removed  it.  The 
error  in  the  estimate  can  be  shown  to  be 


Error  =  LC„L^  , 

and  the  error  propagator  is 


(4-25) 


Error  Propagator  = 


N 


•Trace 


(4-26) 


The  error  propagator  for  the  SWAT  system  is  0.426,  which  means  that  less  than  half  the  phase 
variance  due  to  tilt  jitter  is  effective  in  causing  a  figure  variance  of  the  stitched  beam.  With  these 
results  the  performance  of  a  4-beacon  system  at  zenith  and  at  45  degrees  off  zenith  arc  found  and 
plotted  for  /arious  altitudes  in  Figure  4-9. 
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Figure  4-9.  Predicted  Swat  performance  with  four  synthetic  beacons  as  a  function  of  beacon 
altitude  and  zenith  angle  for  the  HV-21  model.  Stitching  ^ects  have  been  included. 
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5,  SWAT  PERFORMANCE  PREDICTIONS 


The  SWAT  experimental  hardware  has  the  capability  to  probe  an  extensive  range  of  system 
parameters  relevant  to  both  Rayleigh-beacon  and  sodium-beacon  deployment  straiegies.  In  this 
section  the  more  important  of  these  effects  arc  studied  using  the  analytic  and  numerical  simulation 
approaches  previously  described.  In  order  to  permit  these  data  to  be  easily  accessed,  each  of  the 
figures  provided  in  this  section  will  have  a  standard  presentation  format  that  includes  a  chan 
comparing  near- field  variance  for  the  HV-21  model  using  the  two  analysis  procedures  and  a 
second  graph  containing  Strchl  predictions  for  the  HV-21  and  SLC-Day  models  obtained  from  the 
numerical  simulation.  In  aU  of  the  cases  presented,  the  agreement  between  the  two  analysis 
procedures  is  found  to  be  excellent. 

Default  parameters  for  the  next  set  of  performance  curves  are  summarized  in  Table  5-1. 
Unless  otherwise  indicated,  the  pertinent  analysis  inputs  are  as  specified  in  this  list  In  all  cases  it 
should  be  assumed  that  turbulence  effects  dominate  performance . 


TABLE  5-1 

Baseline  Simulation  Parameters 


Atmospheric  Model : 

HV-21 

Aperture  Diameter : 

55  cm,  no  obscuration 

Beacon  Wavelength ; 

500  nm 

HEL  Wavelength ; 

500  nm 

Zenith  Angie : 

OP 

Source  Geometry : 

Single  Source 

Source  Range : 

5  km 

Source  Diameter ; 

Point  Source 
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5.1  SINGLE-BEACON  PERFORMANCE  ESTIMATES 

The  first  series  of  SWAT  experiments  will  study  the  degree  of  turbulence  correction  that  can 
be  achieved  with  a  single  Rayleigh  beacon.  With  this  geometry  the  main  contributors  to  the  figure 
variance  will  be  the  unsampled  atmosphere  above  the  beacon,  focal  anisoplanatism  below  the 
beacon,  and  anisoplanatic  errors  due  to  beacon  misplacement 

5.1.1  Beacon-Altitude  Dependence 

The  upper  chart  in  Figure  5-1  shows  the  expected  figure  variance  as  a  function  of  beacon 
altitude  between  2  and  15  km  for  a  point  source  at  zenith.  The  numerical  and  analytical  results  are 
given  by  the  dashed  and  solid  curves  respectively.  The  agreement  between  the  two  computational 
approaches  is  seen  to  be  extremely  good  over  the  entire  parameter  range. 

The  flguie-Strehl  predictions  derived  from  the  numerical  study  are  summarized  in  the  lower 
chart  of  Figure  5-1,  where  the  HV-21  and  SLC-Day  results  have  been  represented  by  the  dashed 
and  solid  curves  respectively.  Note  that  the  two  models  show  a  cross-over  point  at  the  baseline 
altitude  of  5  km,  and  differ  by  only  30%  for  a  15  km  scattering  altitude.  On  the  basis  of  this 
analysis,  one  would  place  a  rough  upper  limit  of  0.5  on  the  achievable  Strehl  for  single-beacon 
Rayleigh  correction  for  scattering  altitudes  between  5  and  10  km. 

5.1.2  Zenith-Angle  Dependence 

As  discussed  in  Section  4.2,  both  the  upper-altitude  and  lower-altitude  figure-error 
components  associated  with  synthetic  beacons  have  an  explicit  sec(^)  dependence  for  a  fixed 
source  altitude  H.  For  a  fixed  range,  L,  the  error  variance  displays  a  faster  angular  dependence 
due  to  the  cos((^)  reduction  in  beacon  altimde.  These  effects  are  iDustrated  in  Figure  5-2  for  a 
5-km  beacon  range  and  zenith  angles  between  0°  and  60°. 

5.1.3  Beacon-Offset  Dependence 

One  of  the  many  experiments  proposed  in  the  original  SWAT  test  plan  involves  the 
displacement  of  the  beacon  from  its  nominal  point-ahead  position,  as  a  means  of  measuring  tilt 
anisoplanatic  effects.  Although  6^  is  about  10  firad  for  both  the  HV-21  and  SLC-Day  models. 
Figure  5-3  shows  that  essentially  no  degradation  is  predicted  for  beacon  offset  angles  as  large  as 
15  jirad.  This  response  can  be  attributed  to  an  effect  that  might  be  labeled  the  "fractional 
isoplanatic  angle",  which  is  related  to  the  5/3  turbulence  integral  between  the  transmitter  and  the 
beacon  altitude.  For  the  baseline  conditions  in  which  a  5  km  source  is  used  to  probe  HV-21  a 
turbulence  profile,  the  fractional  isoplanatic  angle  is  approximately  25  prad. 
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Figure  5-1.  Single-beacon  performance  as  a  function  of  beacon  altitude.  Calculations  are 

performed  at  zenith  for  a  500  nm  beam. 
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STREHL  ERROR  VARIANCE  (rad 


Figure  5-2.  Single-beacon  performance  as  a  function  of  zenith  angle  for  a  beacon  range  of  5  km. 
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STREHL  ERROR  VARIANCE  (rad 


Figure  5-3.  Single-beacon  performance  as  a  function  of  source  displacement.  The  beacon 

altitude  is  5  km  and  the  zenith  angle  is  0  . 
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5.1.4  Beacon-Size  Dependence 

The  generation  of  synthetic  beacons  using  Rayleigh  or  resonance  backscatter  requires  a  non¬ 
trivial  expenditure  of  laser  energy,  and  this  task  can  be  greatly  simplified  if  beams  of  relatively 
poor  quality  can  be  employed.  Following  the  same  logic  applied  above,  one  would  expect  that 
beam  diameters  of  the  order  of  the  fractional  isoplanatic  angle  would  provide  adequate 
performance.  This  effect  is  demonstrated  in  Figure  5-4,  in  which  virtually  no  change  in  figure 
variance  is  indicated  for  beacon  profiles  as  large  as  60  cm.  This  result,  which  has  been  confirmed 
by  an  independent  study,  predicts  that  compensated  beacons  will  not  be  a  requirement  for  the 
SWAT  experiments.  The  current  design  of  the  SWAT  dye-laser  system  assumes  that  beam  quality 
values  in  the  range  of  10  to  20  are  tolerable. 

5.2  PERFORMANCE  ESTIMATES  FOR  FOUR  BEACONS 

As  discussed  in  Section  2  of  this  report,  the  use  of  multiple  synthetic  beacons  has  been 
proposed  as  a  means  of  reducing  low-altitude  focal  anisoplanatism,  thus  allowing  the  synthetic 
beacon  concept  to  be  applied  to  apertures  of  arbitrary  extent.  The  SWAT  system  is  equipped  to 
deploy  as  many  as  five  dye-laser  beams  in  rapid  sequence.  This  capability  will  permit  direct 
comparisons  of  single  and  multiple-beacon  performance  for  a  variety  of  altitudes,  zenith  angles, 
and  source  displacements.  The  two  most  important  tests  to  be  conducted  with  this  four-beacon 
geometry  are  described  below. 

5.2.1  Beacon-Altitude  Dependence 

As  shown  in  Figure  2-2,  the  four-beacon  SWAT  experiment  will  place  one  source  over  the 
center  of  each  of  four  aperture  quadrants.  In  the  simulation  process,  phase  gradients  are  derived 
from  the  associated  artificial  source  within  each  quadrant,  and  overall  tilt  measurements  of  each 
beacon  are  made  as  a  means  of  removing  anomalous  section  tilts.  Finally,  a  standard  Gauss 
estimator  is  used  to  convert  the  corrected  gradients  into  a  single  phase  array.  The  expected  figure 
variance  and  Strehl  ratio  values  related  to  this  process  are  given  in  Figure  5-5  as  a  function  of 
beacon  altitude. 

The  results  shown  in  this  figure  can  be  directly  compared  to  the  single-source  calculations 
presented  in  Figure  5-1.  For  the  baseline  altitude  of  5  km,  a  Strehl  improvement  of  about  25%  is 
obtained  for  both  the  HV-21  and  SLC-Day  models. 
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Figure  5-4.  Single-beacon  performance  as  a  function  of  beacon  diameter.  The  beacon  altitude  is 

5  km  and  the  zenith  angle  is  0°. 
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STREHL  ERROR  VARIANCE  (rad 


Figure  5-5.  Four-beacon  performance  as  a  function  of  beacon  altitude.  Calculations  are  per¬ 
formed  at  zenith  for  a  500  nm  HEL. 
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5.2.2  Error  Due  to  Beacon  Misplacement 

The  purpose  of  the  section-stitching  algorithm  is  to  render  the  gradient  measurements 
insensitive  to  small  changes  in  relative  beacon  placement  due  to  atmospheric  effects  or 
imperfections  in  the  beam-placement  system.  Figure  5-6  shows  the  performance  as  a  function  of 
the  magnitude  of  the  beacon-placement  error  with  respect  to  the  center  of  the  apenure  section. 
Note  that  essentially  no  change  is  seen  for  placement  errors  as  large  as  ±15  cm. 

5.3  DUAL-ALTITUDE  GEOMETRIES 

One  way  to  obtain  the  benefit  of  a  high-altitude  source  while  retaining  an  acceptable  signal 
level  is  to  perform  the  phase  correction  in  two  distinct  stages:  the  first  stage  entails  a  conventional 
Rayleigh  measurement  using  resized  subapenures  that  provide  high-spatial-frcquency 
information,  and  the  second  entails  section-tilt  measurements  by  large  subapenures  viewing  the 
returns  from  a  more  distant  beacon.  This  so-called  dual-altitude,  or  hybrid  approach  can  be 
implemented  with  illumination  sources  of  relatively  modest  power.  A  0.5-J  laser  pump  at  the 
sodium  resonance  wavelength  has  already  been  constructed  for  this  purpose.^ 

5.3.1  Single  Rayleigh  Beacon  with  One  Sodium  Beacon 

The  physical  structure  of  the  high-altitude  gradient  subapenures  need  not  be  matched  to  the 
number  of  beacons  used  in  the  low  altitude  array.  The  high-altitude  measurements  serve  to  provide 
tilt  corrections  over  contiguous  aperture  sections  that  can  be  arbitrarily  specified.  In  the  analysis 
described  in  Figure  5-7,  section-tilt  measurements  derived  from  a  90-km  source  were  applied  to 
apenure  quadrants  subsequent  to  a  phase  reconstruction  based  on  low-altitude  gradients  obtained 
from  a  single  beacon.  The  results  are  given  as  a  function  of  the  position  of  the  low-altitude  source. 
When  compared  with  Figure  5-1,  it  can  be  seen  that  this  hybrid  scheme  provides  a  45%  Strehl 
improvement  over  the  performance  of  a  single  beacon  at  5  km. 

5.3.2  Four  Rayleigh  Beacons  with  One  Sodium  Beacon 

Dual-altitude  turbulence  sampling  was  originally  proposed  as  a  means  of  removing  the 
residual  section-tilts  associated  with  multiple-beacon  stitching  algorithms.  Figure  5-8  gives 
performance  estimates  for  a  hybrid  scheme  in  which  an  initial  phase  is  obtained  from  a  four-beacon 
Rayleigh  geometry,  followed  by  section-tilt  correction  using  gradient  measurements  provided  by  a 
sodium  source.  This  five-source  compensation  process  yields  roughly  a  20%  higher  Strehl  than 
the  hybrid  approach  in  which  a  single  Rayleigh  beacon  is  deployed. 
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Figure  5-6.  Performance  of  a  four -beacon  system  with  one  beacon  displaced  from  the  center  of 
its  section.  The  displacement  is  radial  with  respect  to  the  center  of  the  full  aperture. 
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1.0 


Figure  5-7.  Analysis  of  a  hybrid-beacon  geometry  in  which  the  high-spatial-frequency  informa¬ 
tion  is  derived  from  a  single  beacon  placed  at  the  specified  altitude  and  tilt  correction  is  subse¬ 
quently  applied  to  aperture  quadrants  based  on  returns  jrom  a  single  90-km  source. 
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Figure  5-8.  Analysis  of  a  hybrid-beacon  geometry  in  which  the  high-spatial-frequency  informa¬ 
tion  is  derived  from  a  set  of  four  beacons  placed  at  the  specified  altitude  and  tilt  correction  is 
subsequently  applied  to  the  four  aperture  sections  based  on  returns  from  a  single  90-km  source. 


5.4  PERFORMANCE  COMPARISONS  FOR  THE  PRINCIPAL  BEACON 
GEOMETRIES 

A  summary  of  the  performance  results  developed  in  this  section  is  presented  in  Figures  5-9 
and  5-10.  Both  illustrations  provide  Strehl  estimates  for  the  HV-21  and  SLC-Day  models  based 
on  output  from  the  numerical  code  ELVIRA.  The  four  principal  beacon  geometries  are  described 
for  a  range  of  altitudes  and  zenith  angles.  These  constructs  include  single-altitude  one  and  four- 
source  arrays,  along  with  the  associated  hybrid  schemes  in  which  a  sodium  source  is  used  to  sense 
residual  section  tilts. 

A  number  of  general  observations  can  be  made  on  the  basis  of  the  curv'es  shown  in  Figures 
5-9  and  5-10.  These  comments  have  direct  relevance  to  experimental  priorities  for  future  SWAT 
field  tests. 

1)  The  Strehl  curves  are  not  strongly  model  dependent  and  probably  represent 
worst-case  examples  for  night-time  operation  at  the  Maui  observatory.  It  should 
be  noted,  however,  that  all  calculations  have  assumed  high  signal-to-noise  m  both 
the  wavefront  sensor  and  far-field  scoring  camera. 

2)  Far-field  performance  at  small  zenith  angles  is  expected  to  be  relatively 
insensitive  to  beacon  altitude  between  7  and  15  kra.  Above  7  km  most  of  the 
strong  turbulence  is  well  characterized  and  focal  anisoplanatic  effects  for  a  60-cm 
aperture  are  not  severe.  For  this  reason  it  is  likely  that  the  5-7  km  range  will  be 
be  of  most  interest  in  order  to  maintain  a  strong  sipal  level. 

3)  An  enhancement  in  Strehl  between  the  four-beacon  and  single-beacon 
configuration  is  likely  to  be  difficult  to  measure.  On  this  point  the  results  from 
the  two  turbulence  models  are  somewhat  different,  but  under  no  conditions  is  the 
improvement  expected  to  be  larger  than  30%.  In  order  to  detect  a  factor  of  this 
magnitude,  the  SWAT  hardware  will  need  to  be  operating  close  to  its  theoretical 
limit,  and  a  large  sample  ensemble  will  be  required  to  provide  good  measurement 
statistics.  There  appears  to  be  a  significantly  higher  probability  of  distinguishing 
between  the  single-altitude  and  hybrid-sensing  geometries. 
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Figure  5-9.  Comparison  of  far -field  performance  for  the  four  major  beacon  geometries  sched¬ 
uled  for  test  by  the  SWAT  experimental  system  as  a  function  of  the  postion  of  the  low-altitude 
beacon  array.  Quadrant-tilt  correction  is  derived  from  a  single  beacon  at  90  km. 
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Figure  5-10.  Comparison  offar-field  performance  for  the  four  major  beacon  geometries  sched¬ 
uled  for  test  by  the  SWAT  experimental  system  as  a  function  of  zenith  angle.  Quadrant  tilt  cor¬ 
rection  is  derived  from  a  single  beacon  at  90  km. 
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In  addition  to  the  Rayleigh  and  hybrid-beacon  designs  covered  in  this  section,  the  SWAT 
program  is  also  funding  the  construction  of  a  7-joule  sodium-frequency  laser  based  on  the  sum- 
frequency  concept^,  which  could  be  used  to  generate  a  single  source  at  90  km.  At  this  altitude, 
focal  anisoplanaiic  effects  become  vinually  negligible  for  a  60-cm  aperture.  Performance  curves 
showing  the  zenith  angle  dependence  for  a  90-km  source  are  given  in  Figure  5-11.  It  can  be  seen 
that  Strehl  ratios  in  excess  of  80%  are  expected  for  angles  as  large  as  45°. 

It  should  be  emphasized  at  this  point  that  the  most  imponant  output  of  the  SWAT  experiments 
will  be  a  validation  of  the  basic  HAVE  REACH  concept  and  a  verification  that  single-beacon  Strehl 
ratios  close  to  the  predicted  values  can  be  consistendy  achieved.  The  exploration  of  a  more  general 
parameter  space  involving  sources  of  varying  size,  configurations,  and  altitudes  will  require 
considerably  more  care  in  order  to  detect  subde  differences  in  performance.  It  will,  however,  be 
important  to  demonstrate  good  correlation  between  experiment  and  theory  for  a  wide  range  of 
synthetic-beacon  constructs  in  order  to  allow  performance  extrapolations  for  apertures  of  larger 
dimensions  to  be  made  with  high  confidence. 
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Figure  5-1 1 .  Far-field  performance  as  a  function  of  zenith  angle  for  a  single  sodium  beacon. 


6.  SUMMARY 


As  part  of  a  comprehensive  analysis  effon  in  suppon  of  the  SWAT  experimental  program, 
independent  numerical  and  analytical  techniques  have  been  developed  to  provide  estimates  of 
system  performance.  With  the  aid  of  this  software  it  is  now  possible  to  accurately  predict  both  the 
near-field  phase  variance  or  far-field  Strehl  for  a  wide  range  of  synthetic  beacon  geometries. 
Despite  the  dissimilarity  of  these  two  computational  approaches,  the  agreement  between  the  figure- 
variance  estimates  is  extremely  good  for  all  of  the  scenarios  studied. 

In  the  system's  baseline  configuration,  in  which  a  single  beacon  is  deployed  at  a  range  of  5 
km,  both  the  SLC-Day  and  HV-21  atmospheric  models  yield  Strehl  predictions  in  excess  of  0.3. 
This  number  exceeds  the  uncorrected  Strehl  by  nearly  a  factor  of  100,  and  the  measurement  of  an 
improvement  of  this  magnitude  would  provide  definitive  proof  of  the  efficacy  of  the  synthetic- 
beacon  concept. 

The  intent  of  the  multiple-beacon  stitching,  hybrid-beacon,  and  sodium-beacon  experiments  is 
to  test  the  scalability  of  these  approaches  to  apertures  of  much  larger  diameter.  Although  a 
measurable  Strehl  improvement  is  of  marginal  probability  for  the  first  two  of  these  constructs,  the 
single-sodium  experiments  are  likely  to  provide  a  clear  verification  of  the  benefits  of  high-altitude 
sampling. 
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APPENDIX  A 

FOCAL  ANISOPLANATISM 


A.l  THEORY 

One  can  use  Mellin  transfonn  techniques  to  calculate  the  phase  variance  of  a  beacon  at  a  given 
altitude,  H,  with  all  the  phase  present  using  Equation  (3.7.4)  of  Technical  Repon  807^  which  is 

<t2=0.5/:oV5/3(^]  = 

Often  one  is  interested  in  the  phase  variance  with  piston  or  piston  and  tilt  removed.  The  phase 
variance  due  to  piston  and  tilt  can  be  calculated  separately  and  subtracted  from  the  total  phase 
variance  to  obtain  these  variances.  The  filter  functions  to  remove  piston  and  tilt  are  given  in 
Equations  (3.4.20)  and  (3.4.23)  and  the  remark  at  the  bottom  of  page  102  of  TR-807.  The  piston 
removal  filter  function  for  a  beacon  at  a  distance  L  is 
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The  filter  function  to  remove  the  phase  variance  due  to  tilt  is 
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These  must  be  inserted  into  the  standard  expression  given  below  to  find  the  variance 

L 

(p-=Q.2Qil3k^\dzCl{z)\dKf(K)F^{K,  z)  . 

0 


(A-2) 


(A-3) 


(A-4) 


To  find  the  piston  variance  insen  Equation  (A-2)  into  (A-4),  perform  the  angular  integration,  and 
obtain 


J^{kDI2)  J^[kD{\-zI  DU) 
kDI2  KD{\-ziL)l2 

Similarly,  the  tilt  variance  is 


cil  =  5.2\2ki]dzCl{z)]dKK~^‘'^ 
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L  »» 

cif  =  lQ.%5ki\dzClU)  \dKK-^''^ 

0  0 

Evaluating  the  integral  in  kappa  space  poses  some  difficulty.  The  two  terms  in  brackets  almost 
cancel  especially  for  very  high  beacon  altitudes.  If  the  bracket  is  expanded,  three  terms  result. 
Two  of  these  terms  can  be  evaluated  by  table  lookup  for  the  tilt  case.  For  piston  each  of  these 
terms  is  divergent.  The  term  that  is  the  product  of  the  two  terms  must  be  evaluated  using  the 
Mellin  convolution  integral.  To  obtain  sufficient  accuracy  a  large  number  of  terms  in  the  resulting 
series  expansion  v^ould  have  to  be  retained  to  get  an  accurate  final  answer  for  the  tilt  case.  A  more 
productive  way  to  manipulate  the  integral  uses  the  multiplication  formula  for  Bessel  functions 
given  in  Equation  (8.535)  from  Gradshteyn  and  Ryzhik*  which  is. 

Il-A^kl.  (A-7) 

For  this  problem 


jk=0^‘ 


J-ijKD/l)  J^jKDiy-zl  DU) 
kDH  kD{\-zIL)I2 


(A-6) 


A=l-r/L  , 


(A-8) 


which  is  less  than  unity  so  that  the  inequality  is  satisfied.  The  first  three  terms  of  the  expansion  arc 

,2' 


yv(Ajc)  =  A* 


Jvix)  +  y  v+i  (X)^(l  -ZI1L)  +  -  2  /  2L)| 


(A-9) 


The  term  in  brackets  to  find  the  piston  in  Equation  (A-5)  becomes 

J2{x){\-zi: 


J\{x)  Ji{x[l-zlL])  z 
X  xil  —  z/L)  L 


AW, 

\xz 

7 

u 

(A- 10) 


The  terras  that  diverge  exactly  cance'  in  this  approach.  The  term  in  brackets  to  find  the  tilt  in 
Equation  (A-6)  is 


J2{x)  J2{x{\-zIL]) 

X  x{\-zlL) 


f,  32  ,  2^  I  Ji{x)X2 


2  L 


(A-11) 


The  entire  first  term  does  not  cancel  in  this  approach.  The  reason  why  there  is  exact  cancellation  of 


70 


the  piston  term  and  not  the  tilt  term  is  that  piston  measured  at  any  altitude  adds  up  with  the  other 
piston  contributions  to  get  the  total  piston.  The  tilt  measured  at  any  altitude  is  magnified  by  the 
ratio  of  the  diameters  when  it  is  added  to  the  other  tilts;  therefore,  there  is  not  perfect  cancellation  of 
the  tilt.  The  first  term  is  by  far  the  most  significant  for  both  piston  and  tilt.  If  these  terms  are  put 
back  into  the  variance  expression  and  squared  as  required,  the  Mellin  transforms  can  be  performed 
to  give  the  following  approximation  for  the  variances. 


The  piston  variance  is 

The  tilt  variance  is 


0.0835^+0.01378^  +  0.00697^-0.0591^  +  0.045^ 
l}  1}  U  . 


(A- 1 2) 


0.355^-0.484^  +  1.00^-0.543^  +  0.116^ 
Is  IS  IS 


(A-13) 


The  piston  removed  variance  is 


0.348D 


x5/3 


• 


(A-14) 


The  piston  and  tilt  removed  variance  is 


= 


^0.348D^^^^ 


V  J 


-<^p  -of  • 


(A-15) 


A-2  SUMMARY 

A  simple  way  to  compute  the  phase  variance  with  focal  anisoplanatism  with  and  without  the 
presence  of  piston  and  tilt  has  been  provided.  This  is  easily  extended  to  remove  other  Zemike 
components.  A  Basic  computer  program  has  been  written  for  the  Macintosh  that  performs  these 
calculations. 


Using  the  above  results,  a  good  approximation  to  the  piston  removed  phase  is 


/^5/3 


0.17^-0.02756^-0.0139^  +  0.118^-0.09 


U  L  Lf 

Similarly,  a  good  approximation  to  the  piston  and  tilt  removed  phase  is 


(A- 16) 
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^-0.877^  +  0.94%-2.014^  + 1.204^-0.322^ 
Is  ^  U"  U  L7  u  L° 


(A- 17) 
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APPENDIX  B 

FOCAL  ANISOPLANATIC  TILT 


B-l  INTRODUCTION 

Consider  two  beacon  sources  at  an  altitude  L  with  a  separation  b  as  shown  in  Figure  Bl. 
Let  the  position  of  each  of  these  beacons  be  measured  with  a  common  aperture.  The  tilt  difference 
measured  between  these  two  beacons  will  be  equal  to  the  real  difference  between  these  two 
beacons  and,  in  addition,  there  will  be  a  stochastic  difference  with  zero  mean  that  is  due  to  the 
different  paths  that  the  two  beacon  rays  take  through  the  atmospheric  turbulence.  It  is  the  two-axis 
rms  stochastic  tilt  difference  between  these  two  point  beacons  that  is  calculated  in  this 
memorandum. 


Figure  Bl .  Focal  anisoplanatic  geometry. 


B.2  THEORY 

The  tilt  difference  between  two  infinite  plane  waves  looking  in  different  directions  with  infinite 
outer  scale  size  has  been  treated  in  Section  3.6.2  Md  3.6.3  of  TR-807.  This  effect  is  referred  to  as 
tilt  anisoplanatism.  In  this  memorandum,  rather  than  having  infinite  plane  waves,  one  is  dealing 
with  focused  waves.  The  only  difference  between  this  problem  and  the  one  that  has  been  solved  is 
that  the  tilt  filter  function  must  be  the  one  appropriate  for  focused  beams.  The  filter  function  is  a 
product  of  the  term  for  artisoplanatism  and  the  one  to  extract  tilt  of  a  beam  focused  at  L.  The  mean 
squared  tilt  variance  can  be  written  down  as 


=  0.2073/^1  dz  Cliz)  J  die  k 


r;^-ll/3 


16  f  f  J2{icD[l-z/L]/2y 


k^D 


kD[1-z/L]/2 


1  -cos|>c-  ^(z)| 


(B-l) 


where 
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The  angular  integration  can  be  performed  and  if  the  substitutions 

(B-3) 
(B-4) 

arc  made,  then  the  tilt  variance  is 


and 


x=  /cD[l-z/L]/2  , 


2b  z 


rl  =  dzcli.2){\ - z / - Joix / ;)]  .  (B-5) 

^0  0 

In  the  last  formula  the  zenith  dependence  was  included.  It  was  assumed  that  as  the  zenith  anglp 
changes  the  distance  of  the  beacon  to  the  transmitter  changed  to  keep  the  beacon  altitude  constant. 
To  get  the  zenith  dependence  with  that  assumption  the  substitutions  2-4zsec(<^)  and 
L  — >  Lsec(^)  were  made.  For  large  values  of  t  the  last  integral  goes  to  zero.  This  occurs  at  low 
altitude  because  there  is  very  little  separation  between  the  rays  of  the  two  beacons.  There  is  a  large 
separation  between  the  rays  close  to  the  beacon  but  there  is  also  very  little  contribution  from  this 
region  because  of  the  demagnification  factor  of  the  tilt.  This  is  seen  from  the  term  on  the  right  of 
the  turbulence  strength.  The  Mellin  convolution  integral  can  be  used  to  conven  the  last  integral  into 


/ 


_l _ 1_ 

2^fn  2m 


jdsi2t)~^^r 

c 


S  +  ^,  — 5  +  y,  —5* 
-5  +  ^,-5  +  ^,  J  +  1 


(B-6) 


Since  Zi  =  0,  the  path  of  integration  should  be  closed  in  the  left-half  plane  for  small  values  of  the 
parameter  and  in  the  right  half  plane  for  large  values.  The  parameter  in  the  integral  is  large  for  low 
altitude  and  small  for  high  altitudes;  therefore,  to  get  a  rapidly  converging  series  one  wants  one 
power  series  for  high  altitudes  and  the  other  power  series  for  low  altitudes.  The  power  series  has 
poles  at  5  =  0,  and  s  =  -1/6  -  n  for  n  =  0,  1,  2,  ...  .  The  parameter  is  equal  to  unity  at  an 
altitude  of 


L 

b/D-^-l  ' 


(B-7) 


From  the  pole  locations,  the  power  series  solution  for  high  altitudes  can  be  written  as 


/  = 


0.2052- 


n=0 


n! 


n  +  4,  «  +  -5' 
n  +  5,  «  +  3,  -  n  +  T 


(B-8) 


The  power  series  good  for  low  altitudes  is 

2V^/  = 


-£(=fc,)-r 

«=1 


«  +  i,-n  +  7 


£lzl):(2rr«^^V 


n=0 


«+■!, 


n-3 


„  ,  5  „  ,  1  „  ,  10 


(B-9) 


The  transition  altitude  when  one  should  change  from  one  power  series  to  the  other  is  r, ,  however, 
for  good  convergence  with  only  a  few  terms  use  the  second  power  series  whenever  r  >  z,  /  2  and 
the  first  power  series  for  lower  altitudes.  Substitute  these  results  back  into  Equation  (B-5).  The 
first  few  terms  of  each  solution  give  the  till  variance  as 


14/3 


(1  -  z  /  L)-^  +  0-339^1^  J(1  -  z  / 


12.16 
■^“TiTT 


\dzCl{z) 


z,l2 


vl/3 


(1  - z / - 0.665?!^^^  {\-z/LY- 0.003081 


(fr<- 


iLY 


-3.06xl0'^f^J  (l-z/L)^-6.15xlO"^(-^j  (1-z/L)* 


vl9/3 


(B-10) 


Therefore,  to  evaluate  this  expression  for  a  particular  beacon  separation  and  height,  apenure 
diameter,  and  turbulence  model,  one  has  to  perform  the  numerical  integration  over  the  turbulence 
profile.  A  Basic  program  has  been  written  to  perform  this  calculation.  As  an  example,  for  a 
beacon  altitude  of  100  kra  and  separation  of  0.3  m  with  a  0.6-m  receiving  aperture  the  tilt  standard 
deviation  is  42,  41,  and  72  nanoradians  for  the  SLC-Day,  SLC-Night  and  HV-21  turbulence 
models. 
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B-3  SUMMARY 

Mellin  transform  techniques  have  been  used  to  denve  an  expression  for  the  vanance  of  the  ult 
difference  between  two  beacons  that  emanate  from  point  sources  at  the  same  altitude  and  whose 
positions  are  measured  by  a  common  apenure.  A  computer  program  is  available  to  perform  this 
calculation. 
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APPENDIX  C 

CORRELATION  FUNCTIONS  OF  TILT  AND  OF  THE  TILT 
AND  PISTON  COMPONENTS  OF  FOCAL  ANISOPLANATISM 


C-1  INTRODUCTION 

Consider  the  tilt  or  piston  difference  between  a  collimated  beam  and  a  focused  beam.  In  this 
memorandum  the  correlation  function  for  these  tilt  and  piston  components  is  calculated  first  for  the 
case  in  which  the  point  source  is  directly  over  the  apenure  and  then  for  the  case  in  wliich  the  point 
source  is  displaced  from  the  central  location.  The  correlation  function  for  tilt  is  a  function  of  the 
direction  of  the  displacement  relative  to  the  tilt  component.  Surprisingly,  it  is  found  that 
perpendicular  components  of  tilt,  which  are  uncorrelated  in  a  given  aperture,  are  correlated  for 
aoertures  that  are  displaced  in  a  skewed  direction. 

In  focal  anisoplanadsm  the  major  part  of  the  tilt  component  comes  about  because  the  tilt 
measured  at  higher  altitudes  is  sampled  over  a  smaller  area  than  the  collimated  beam  and  is 
demagnified  as  the  beam  expands.  The  result  is  that  the  tilt  and  the  tilt  component  of  focal 
anisoplanadsm  to  first  order  have  the  same  spadal  statistics  and  subsequently  the  same  correlation 
function. 


C.2  CORRELATION  FUNCTION  FOR  A  CENTERED  POINT  SOURCE 

For  some  problems  it  is  necessary  to  compute  the  correlation  function  of  the  tilt  and  piston 
variance  components  of  focal  anisoplanadsm.  Using  the  filter  function  method  it  is  relatively  easy 
to  calculate  this  correlation.  The  correlation  function  for  the  piston  and  tilt  is  defined  to  be 


2 


=  !-■ 


(C-I) 


The  denominators  of  the  equations  have  been  calculated  in  the  appendix  on  focal  anisoplanadsm. 
The  numerators  can  be  calculated  from  the  expressions  in  that  memorandum  by  incorporating  the 
anisoplanatic  filter  function.  There  are  three  correlation  functions  of  the  tilt  that  must  be  calculated. 
The  numerator  terms  can  be  written  as 


77 


L 

=  \.659l4\dzCl{z)\dKK~^^'^ 

0 


j^jKon) 

kDJ2 


Jy{KD{\-zl  L)ll) 
kD{\~zIL)I1 


1-cos 


(f  5) 


(C-2) 


and  for  the  a  and  h  components  of  tilt  it  is 


^lab-  =  +  d)f  =  2Mki\dzCl{z)  J dK \g,^  -  G,^  exp[yic-  d|  .  (C-3) 


where  the  complex  filter  function  applies  to  that  component  of  tilt.  The  correlation  coefficient 
between  the  jc-components  of  dlt  is 


2 

^^ca'=[crrr(n-<^a(^  +  ^)]  =  6.636^  JizC^Cz)  Jdic  V 


11/3 


j^Xkdh)  j^jKDG-zmfi) 
kDH  ”  kD{\-zIL)I2 


-a 


cos^(0)  l-cos(/c  d) 


(C-4) 


The  correlation  function  of  the  x-tilt  is  the  correlation  of  the  x-tilt  component  with  the  x-tili 
displaced  in  an  arbitrary  direction.  The  formula  for  the  y  component  will  be  written  by  inspection 
after  this  derivation.  A  useful  integral  in  the  evaluation  is 

In  n 

f  d6  cos( 6)  sin( 6) cos[ >cd  cos( 0)]  =  2j d6 cos( 6)  sin(  0) cos[  Kd  cos{  0)]  =  0  .  (C-5 ) 

0  0 

If  the  last  integral  is  broken  up  into  a  part  from  0  to  7t/2  and  from  /t/2  to  n,  and  the  substitution 
$=  9-  nis  made  in  the  second,  the  two  pans  cancel.  If  the  angle  between  the  displacement 
and  the  a:-axis  is  0',  the  angular  integration  can  be  expanded  using  the  trigonometric  identity 

cos(fl  +  0)  =  cos(a)cos(f?)-sin(a)sin(f7)  ,  (C-6) 


to  give 
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2n 


j4i0cos^(0)|l-cos[K'cfcos(0  +  ^■)||  =  jd0COS^i<P  -  0,)|l  -cos[K'i/cos(0)]| 

0  0 

=  Jii0|cos^(0,)cos^(0)  +  sin^(fl,)sin^(0)-2cos(0i)cos{0)sin(0j)sin(0)j 
0 

x|l-cos[vdcos(0)]J.  iC-l) 

The  integral  of  the  last  term  in  brackets  is  zero  since  it  is  the  same  integral  as  in  Equation  (C-5). 
The  angular  integration  of  the  second  term  can  be  performed  using  Equation  (3.915  #5)  from 
Gradshteyn  and  Ryzhik  which  is 

Jexp[i^cos(jc)]sin^‘'(x)dx  =  r[v4-^] 7^,(J3);  Rev>-4  (C-8) 

Use  the  trigonometric  identity 


cos^{<p)  =  l-sin^((j3)  (C-9) 

to  express  the  first  term  in  the  form  of  Equation  (C-?)  which  integrals  can  then  be  performed  to 
give 


N, 


Da 


J2{kD/2) 

kDH 


J2{kD{\-zI  L)ll) 
x-Z?(l-z/L)/2 


[sin^(0,)[l ll-J^{Kd)l  Kd]  +  cos^(0.){[l -JoiKd)] -[\/l-J^{Kd)/  ictf]}] . 


(C-10) 


The  Erst  term  of  the  expansions  of  the  Bessel  function  contained  in  the  brackets  in  Equations 
(C-2)  and  (C-10)  is  the  dominant  term.  With  this  single  term  approximation,  which  is  evaluated 
using  Equations  (A-10)  and  (A-1 1),  the  numerator  terms  can  be  written  as 

L  2  “ 

Np  =  \{)Alki\dzCl{z)^\dKK-^''^jl{KDI2)[\~jQ{Kd)]  ,  (C-11) 

0  ^0 

and  for  the  tilt,  the  axial  integration  can  be  performed  to  give 
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(C-I2) 


where 


k~^'^jI{kDI2) 


\-jQi,Kd) 

1  /  2  -  7j  { Kd)  /  Kd 


(C-13) 


The  y-tilt  correlation  coefficient  is  found  by  rotating  the  above  result  by  90  degrees  to  get 

5 /3 

=  ^[sin^{ei--nl2)h^co%^{6-,-Kll){lT-h]]  .  (C-14) 

For  pure  tilt  the  expression  above  has  the  second  moment  divided  by  the  source  distance  squared 
replaced  by  the  zeroth  moment.  When  the  correlation  function  is  calculated,  this  term  is  cancelled 
by  the  denominator  term;  therefore,  in  this  approximation  where  only  the  first  term  of  the 
expansion  of  the  difference  of  Bessel  functions  is  used,  the  correlation  functions  of  tilt  focal 
anisoplanatism  and  tilt  are  the  same. 


The  correlation  coefficient  of  the  perpendicular  components  of  tilt  will  now  be  found.  Use 
Equation  (C-3)  to  get 


L 

=  3.3I8/tJ  JdzC^  (z)  j  dK 
0 


J2{kDI2) 

KDf2 


J2{KD{\-zlL)f2) 

kD{\-zIL)/2 


cos(0)-sin(0)exp|yic- 


(C-15) 


This  can  be  simplified  to  give 

L 

=  3. 3 1 8*^  Jdz  (z)  J  die 

0 

2 

1  -  2sin(0)cos(0)cos^x’ ■  d  j 


J2{kDI2)  J2{kD{1-zI  L)I2) 
kDI2  kD{\-zIL)I2 


(C-16) 
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The  integral  of  the  last  term  over  angle  is 


In 

/=  Jj0sin(0)cos(0)cos[Ar«icos((^- (9')j  .  (C-17) 

0 

Let  0  — 6'=  d,  then  using  the  trigonometric  identities  for  the  sum  of  angles  one  obtains 
2n 

/  =  Jd0sin(0+  0')cos(0+  0’)cos[>cdcos(0)] 

0 

2n 

=  jd9 cos[ Kd cos(0)]|sin(0' ) cos( 0'  )|l  -  2 sin" ( 0)j  +  sin( 0) cos( 0)[cos“  ( )  -  sin " ( 0’ )  | . 

0 

(C-18) 

Using  Equation  (C-5)  one  sees  that  the  last  term  in  braces  integrates  to  zero  giving 

n 

/  =  2sin(0')cos(0')Jd0[l-2sin^(9)]cos[vdcos(0)]  .  (C-19) 

0 

Using  Equation  (C-7)  twice  and  inserting  this  into  Equation  (C-15)  one  obtains 


'J2{kD/2) 

J.iKDil-z!  DU) 

kD/2 

KD{\-zlDn 

X |l  - 2 sin ( 0' ) cos( 0' )[yo {Kd)-2Ji[Kd)l  Kd^. 


(C-20) 


The  correlation  term  which  is  the  last  term  in  braces  is  not  always  zero.  It  is  zero  if  the 
displacement  is  along  the  x  or  y  axis  or  if  the  wavenumber  is  zero  or  infinitely  large.  Again, 
using  only  the  first  terra  of  the  expansion  for  the  expression  in  brackets  one  obtains 

/;3\5/3 

=  20.85^2 [O.O579  +  sin(20’)(2/,  .  (C-21) 


C-3  EVALUATION  OF  CORRELATION  COEFFICIENTS 

The  expressions  developed  in  the  last  section  will  now  be  evaluated  using  Mellin  traniform 
techniques.  Let  a  =  K‘D/2  and  f3  =  D  /  2d,  Equation  (C*l  1)  and  (C-13)  become 
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hi  Ida  -11/3,2/  J  l--/o(o://5) 

.  =  — «  '  JUa) 

LAJ  ll/2-J^(a/p)P/a\ 

Using  the  Mellin  convolution  integral,  the  integrals  over  a  can  be  written  as 


P  9,/^  -y-rrii  I  /V 


2V^  27ii  I  V 


For  the  first  component  of  tilt 


For  the  second  component  of  liit 


‘  A^dlnii  {d 


S  +  -|,  -r  +  y,  —  5* 
■s  +  ^,-s  +  ^,s  +  \ 


5  +  -^.  -5  +  ^,  -5* 

-S+f  .-5  +  J^,5+l 


-J  +  y,  -r +  y,  5  +  4 


(C-23) 


(C-:4) 


(C-25) 


(C-26) 


Since  id  =  0,  the  direction  of  path  closure  is  determined  by  the  size  of  the  parameter  in  the  integral. 
That  parameter  can  be  less  than  or  greater  than  unity  which  requires  the  solution  for  both  ways  of 
closing  the  integration  path.  If  d  >  D,  the  path  of  integration  is  closed  in  the  left-half  plane,  and 
there  are  poles  at  r  =  0  and  at  s  =  -7/6  -  n  forn  =  0,  1,  2  ...  for  the  piston  integral,  and  at  5  = 
0  and  at  5  =  -1/6  -  n  for  n  =  0,  1,2  ...  for  the  first  tilt  integral  and  at  5  =  1/2  and  at5=-!/3- 
n  for  n  =  0,  1,2  ...  for  the  second  till  integral.  \id<D,  the  path  of  integration  is  closed  in  the 
right-half  plane  and  there  arc  poles  at  5  =  «  for  n  =  1,  2,  3  ...  and  at  5  =  4/3  +  n  for  n  =  0,  1,2 

...  for  the  piston  integral  and  at  5  =  «  for  n  =  1,  2,  3  ...  and  at  5  =  7/3  +  n  for  n  =  0.  1,2  ... 

for  the  first  tilt  integral  and  at  5  =  1/2  +  «  for  «  =  1,  2,  3  ...  and  at  j  =  1 7/6  +  for «  =  0,  1,2 

...  for  the  second  tilt  integral.  Therefore,  the  piston  integral  is  equal  to 
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Using  this  in  Equation  (C-22),  one  finds 
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The  piston  variance  is  found  by  using  the  first  term  of  Equation  (A-lO)  in  (A-5)  to  give 

The  piston  correlation  function  can  then  be  written  as 
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(C-33) 


The  piston  is  the  difference  between  the  piston  of  the  parallel  and  focused  beams.  Both  beams 
sample  the  piston  correctly  and  one  expects  that  there  will  be  very  little  correlation  when  the 
apertures  move  far  enough  apan  that  they  are  no  longer  touching.  The  plot  in  Figure  C- 1  indicates 
that  this  is  the  case.  The  correlation  function  essentially  drops  to  zero  at  that  separation. 


Figure  C- 1 .  Correlation  function  for  the  piston  component  of  focal  anisoplanaiism. 
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For  the  first  component  of  tilt,  the  integral  is  equal  to 
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d<D.  (C-35) 


For  the  second  component  of  tilt,  the  integral  is  equal  to 
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The  single  axis  tilt  variance  can  be  found  by  using  the  first  term  of  Equation  (A-1 1)  in  fA-6)  and 
dividing  by  two,  to  obtain 


(C-38) 


The  tilt  correlation  function  can  then  be  written  as 
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(C-40) 


The  correlation  function  is  independent  of  the  turbulence  distribution.  The  tilt  correlation 
functions  for  parallel,  perpendicular  and  45  degree  aisplacement  are  plotted  in  Figure  C-2. 
Observe  that  the  correlation  function  is  significant  even  for  large  separations.  Using  the  relations 
in  Equations  (C-21)  and  (C-28)  in  (C-1),  one  obtains 

Ct^.(d)  =  17.28sin(20,)(/7--2/,)  .  (C-41) 

If  the  values  of  /j-  and  /j  are  insened  into  the  equation  above,  the  correlation  functions  of  the 
perpendicular  components  of  tilt  are  equal  to 
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Figure  C-2.  Correlation  function  for  tilt. 
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(C-43) 


If  the  displacement  is  parallel  to  one  of  the  tilts  the  correlation  is  zero.  For  a  displacement  of  45 
degrees  to  the  tilts  the  correlation  is  plotted  in  Figure  C-3.  There  is  not  as  much  correlation  as  that 
between  parallel  components  of  tilt. 
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Figurs  C-3.  Correlation  function  of  perpendicular  components  of  tilt  with  the  displacement 

45  degrees  to  the  tilt. 
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APPENDIX  D 

FAR-FIELD  IMAGE  ANALYSIS 


The  processing  methods  used  to  extract  estimates  of  Strehl  from  far-field  imagery  are  typically 
reviewed  at  the  beginning  of  each  new  field  exercise,  and  the  SWAT  experiments  have  been  no 
exception.  Over  the  years  a  wide  variety  of  approaches  have  been  developed  to  deal  with  this 
problem,  all  of  which  display  a  fairly  high  susceptibility  to  data  normalization  enx)rs.  Folio., ing  a 
careful  preliminary  analysis  of  the  data  obtained  during  the  SWAT  laboratory  exercise,  a  new 
approach  has  been  devised  that  appears  to  be  somewhat  better  suited  to  this  task.  The  process 
comprises  two  distinct  stages:  during  the  first  step  the  raw  data  are  carefully  corrected  to  remove 
fixed-pattern  noise  and  stray  light  effects;  in  the  second  step  an  estimate  of  the  peak  irradiance 
value  is  obtained  from  an  energy  distribution  function.  The  essential  elements  of  this  process  are 
outlined  below. 

The  high-resolution  CCD  cameras  used  in  the  SWAT  program  can  provide  two-dimensional 
beam  images  of  extremely  high  quality.  The  raw  output  data  Z(x,y)  can  be  modeled  as  follows 

y)=  y)E{x,  y)  +  Zh{x,  y)  ,  (D-1) 

where  E  and  Z  represent  the  input  optical  itradiance  and  the  electrical  outpui  signals  respectively, 
is  the  camera  dwell  time,  is  the  optical  throughput,  d  is  the  pixel  dimension,  g  is  the  the 
camera  gain,  and  Z^  is  the  fixed-pattern  background.  The  background  is  normally  assumed  to  be 
spatially  varying  but  temporally  static  so  that  the  Zf,  term  can  be  eliminated  by  subtracting  camera 
frames  taken  with  the  input  light  blocked.  It  is  further  assumed  that  the  gain  is  independent  of 
pixel  location  but  that  it,  and  the  system  throughput,  cannot  be  accurately  measured.  Therefore, 
the  optical  irradiance  function  at  the  camera  focal  plane  bears  the  following  relationship  to  the 
output  data 

E{x,  y)  =  (1  /  r)[Z{x,  y)-  Z^{x,  y)]  (D-2) 

where  y  is  indeterminate.  Therefore,  a  direct  radiometric  determination  of  the  Strehl  ratio  is 
precluded  but  a  reasonably  accurate  estimate  can  be  obtained  if  the  diffraction-limited  beam  profile 
for  the  collection  aperture  can  be  computed.  In  general,  a  parameter  having  the  form 

Zix,y)-Zf,ix,y) 

2^]^[Z(x,  y)-Zft(x,y)] 

*  y 

will  be  generated  at  some  point  in  the  process,  which  eliminates  the  y  term.  This  energy 
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normalization  is  less  straightforward  than  it  appears,  owing  to  the  diTiculty  of  simultaneously 
achieving  high  spatial  resolution  across  the  beam  while  obtaining  an  accurate  rixasure  of  the  total 
flux.  Separate  narrow  and  wide  field-of*vicw  detectors  can  be  used  for  this  purpose,  but  care  must 
be  taken  to  establish  an  accurate  mutual  calibration  of  the  two  sensors. 

D-l  BASELINE  CORRECTION 

Failure  to  properly  subtract  the  fixed-pattern  background  noise  is  the  greatest  potential  source 
of  error  in  the  Strehl  calculation,  since  a  baseline  offset  on  the  order  of  0.1%  of  the  peak  pixel 
intensity  can  translate  into  a  10%  error  in  the  estimate  of  the  total  beam  energy.  An  offset  error  will 
occur  when  the  ambient  lighting  changes  between  the  data  and  background  frames  or  when  the 
camera  is  exposed  to  signal  light  scattered  from  optical  elements  within  its  field  of  view. 

Recall  that  the  Z-Z^,  operation  is  intended  to  produce  a  data  field  that  is  free  from  fixed- 
pattern  errors.  Assuming  that  the  focused  beam  subtends  only  a  small  fraction  of  the  total  detector 
array  area,  a  histogram  of  Z-Z^  should  display  a  Gaussian  peak  centered  at  zero.  Any 
displacement  of  this  peak  from  the  null  position,  as  illustrated  in  Figure  D-l,  is  evidence  of  a 
baseline  offset  that  should  be  corrected. 

Baseline  compensation  entails  the  formation  of  the  quantity 

Iix,y)  =  Zix,y)  -Zij{x,y)  -Zo  , 

where  Z^  is  the  computed  mean  of  the  baseline  peak  of  the  Z  -  Z^  irradiance  histogram.  While 
the  absolute  value  of  Z^  is  typically  small,  the  Strehl  measurement  error  is  proportional  to  the 
product  of  Zq,  and  the  number  of  pixels  in  the  detector  array.  The  implementation  of  this  simple 
procedure  has  significantly  improved  the  frame-to-frame  uniformity  of  the  Strehl  computations 
performed  on  the  SWAT  laboratory  data. 
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FREQUENCY  of  OCCURRENCE 


Figure  D-1.  Typical  histogram  plot  of  uncorrected  irradiance  at  null. 
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D-2  PEAK  IRRADIANCE  ESTIMATION 

For  a  uniformly  illuminated  circular  aperture,  the  diffraction-limited  beam  profile  is  specified 
by  the  airy  iiradiance  function 


({kDaf)^\ 

j.ikDran'f 

1  j 

kDr/2f 

(D-3) 


where  D  is  the  aperture  diameter,/ is  the  focal  length  of  the  telescope,  and  Q  is  the  pulse  energy - 
It  will  be  convenient  to  approximate  this  function  by  the  Gaussian 
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defines  the  Gaussian  beam  width. 


The  beam  profile  of  an  aberrated  beam  might  be  modeled  in  the  foUowing  way 
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where  S  is  the  peak  iiradiance  (the  usual  definition  of  Strehl)  and  a  is  the  beam  diameter  (which 
will  usually  be  larger  than  cr^).  The  inadiance  can  be  treated  as  a  probability  density  from  which  a 
distribution  function  can  be  formed 
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The  function  F(£)  represents  the  total  focal  plane  energy  attributed  to  those  pixels  exceeding  the 
iiradiance  threshold  E.  F(E)  is  seen  to  be  a  linear  function  of  the  iiradiance,  with  a  slope  that  is 
proportional  to  the  beam  radius  a.  The  x-axis  intercept  Ep  =  QSllnc^  provides  a  means  to 
measure  of  the  peak  iiradiance,  while  the  y-axis  intercept  gives  the  total  energy  included  in  the 
Gaussian  core.  Therefore,  a  linear  fit  to  the  upper  portion  of  the  iiradiance  distribution  function 
win  generate  estimates  of  both  the  Strchl  and  the  effective  beam  diameter.  The  general  appearance 
of  the  F(F)  curves  for  diffraction  limited  and  aberrated  beams  are  illustrated  in  Figure  D-2.  The 
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knee  of  the  curve  will  generally  occur  well  below  the  E^l  point,  so  it  is  usually  safe  to  include  all 
data  points  above  this  value  in  the  linear  regression  operation. 

The  preceding  derivation  provides  a  working  definition  of  the  Strehi  ratio  that  includes  the 
following  analysis  steps: 

1 .  the  formation  of  the  corrected  camera-data  histogram  hil), 

2.  the  creation  of  the  camera  distribution  function  Fil)~^ih(i), 

3.  the  development  of  a  linear  fit  to  the  upper  half  of  the  F(/)  curve  to  obtain  F'(I),  and 

4.  the  determination  of  the  x-axis  intercept  Ip  for  which  F’[}p  j  =  0. 


The  Strehi  ratio  is  now  defined  to  be 


S  = 


2nailp 
*  y 


(D-8) 


As  usual,  the  determination  of  Strehi  requires  an  accurate  measurement  of  the  total  pulse  energy 
and  knowledge  of  the  caiiwra’s  magnification.  Observe  that 
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4  fA 
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where  d  is  the  pixel  dimension,  n  is  the  number  of  phase  sensor  subapenures,  and  d  is  a 
magnification  factor  given  in  units  of  beam  motion  (pixels)  per  wave  of  tilt  per  subaperture.  Since 

A 

the  Strehi  is  proportional  to  d  ,  it  is  important  that  this  magnification  factor  is  remeasured 
following  each  change  to  the  optical  arrangement. 
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Figure  D~2.  Energy  distribution  functions  for  diffraction-limited  and  aberrated  beams. 
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APPENDIX  E 

TURBULENCE  MODEL  PROHLES 


The  Hufnagcl-Vallcy^‘^0  and  SLC-Day  turbulence  profiles  arc  used  in  the  performance  studies 
described  in  this  report  as  a  means  of  indicating  a  range  of  potential  results.  The  Hufnagei-Valley 
profile  actually  refers  to  a  family  of  models  that  have  the  following  general  form 

ChH)  =  5.94xi0~^^(w  /  27)^//^°exp(-//  / 1000)  +  2.7xl0~‘'^exp(-f/  / 1500) 

fE-ll 

+  4exp(-///100) 

where  H  is  measured  in  meters  above  sea  level  and  w  is  the  rms  wind  speed  averaged  over  the  5 
to  20  km  interval.  The  coherence  diameter  is  affected  principally  by  the  low-altitude  constant  A 
(which  is  normally  set  to  1.7x10'*'^  m*^),  while  the  isoplanatic  angle  is  most  sensitive  to  w.  The 
HV-21  model  referenced  to  in  this  report  uses  a  value  of  21  m/scc  for  w;  since  this  results  in  a 
coherence  diameter  of  5  cm  and  an  isoplanatic  angle  of  7  ixrad,  this  version  of  the  model  is  also 
referred  to  as  HV5/7. 


The  SLC-Day  model  is  based  on  the  day-time  Miller-Zieske  profile,  which  can  be 
approximated  by  the  set  of  power-law  segments  given  below 


ciiH)  =  4.01x10'^^//“^°^ 
=  1.3x10"^^ 

=  6.35x10*'' 

=  6.2lxl0~^^ 

This  model  is  also  referred  to  as  SLCSAT  Day. 


19S//<230 
230  <  //  <  850 
850  <  //  ^  7000 

7000  <H<  20000  .  (E-2) 
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